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ABSTRACT 


The  use  of  Doubly-Fed  Induction  Generators  (DFIG)  for  wind  energy  conversion  is 
addressed  in  this  thesis.  It  is  well  known  that  when  the  stator  is  connected  to  the  electric 
grid,  the  rotor  voltage  can  control  both  mechanical  torque  and  reactive  electric  power. 

To  guarantee  efficient  wind  energy  conversion,  it  is  important  to  research  and 
design  more  advanced  control  schemes.  In  this  thesis,  we  first  review  the  basic  theory 
behind  DFIGs  and  Adaptive  Control.  Next  we  design  an  adaptive  controller  for  a  wind 
turbine  using  a  DFIG  and  model  and  simulate  the  system.  In  order  to  create  a  valid 
assessment  on  the  results  of  this  method,  we  compare  the  system’s  perfonnance  with  a 
standard  control  scheme  based  on  proportional  integral  (PI)  controllers  as  proposed  in 
standard  approaches. 
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EXECUTIVE  SUMMARY 


The  possibility  of  fossil  fuel  shortages  in  the  near  future  and  environmental  concerns 
about  power  generation  has  lead  to  an  increasing  demand  for  renewable  energy  (RE) 
installations  worldwide.  Among  them,  wind  energy  (WE)  is  one  of  the  fastest  growing. 
This  fonn  of  energy  production  is  suitable  for  electricity  production  since  wind  is 
available  everywhere.  However,  the  random  nature  of  the  wind  and  the  need  to  provide  a 
stable  source  for  the  grid  calls  for  the  design  of  effective  control  systems  for  wind  energy 
conversion  systems  (WECS).  This  would  improve  both  the  reliability  of  the  energy 
captured  and  the  overall  efficiency  of  the  system. 

Many  large  wind  power  generators  today  are  based  on  the  doubly-fed  induction 
generator  (DFIG).  The  DFIG  is  capable  of  controlling  both  the  active  and  reactive  power 
simultaneously  and  keeping  the  generator  at  its  maximum  efficiency,  even  in  the 
presence  of  wind  speed  variations.  To  take  advantage  of  these  features,  the  DFIG  requires 
proper  control  in  terms  of  its  power  electronics  and  computer  algorithms. 

Standard  controllers,  such  as  proportional  integral  (PI)  controllers  for  pulse  width 
modulation  (PWM)  converters  at  the  rotor  and  source  side  of  the  system  are  the  most 
commonly  used  to  provide  independent  control  of  active  and  reactive  power  output. 

A  different  control  approach  based  on  adaptive  control  theory  is  investigated  in 
this  thesis.  The  motivation  for  this  research  is  the  fact  that  WECS  operate  in 
environments  with  extreme  variations  in  the  operating  conditions,  and  their  operation 
depends  on  the  wind  speed,  which,  in  general  is  a  random  process  that  cannot  be 
controlled.  Adaptive  control  systems,  on  the  other  hand,  have  the  ability  to  “adapt”  to 
variations  in  the  system’s  parameters,  which  is  a  very  attractive  feature.  Among  the 
various  kinds  of  adaptive  control  methods,  the  one  examined  in  this  thesis  is  direct  model 
reference  adaptive  control  (DMRAC).  The  basic  idea  behind  this  approach  is  to  make  the 
system  behave  as  closely  as  possible  to  a  reference  model  in  the  presence  of  uncertain 
and/or  time-varying  dynamics  of  the  system  to  be  controlled.  To  achieve  this,  the  system 
is  modeled  in  tenns  of  a  vector  of  parameters  0c(t)  that  may  be  known  or  approximately 
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known.  This  parameter  vector  is  updated  by  an  adaptive  law  using  input  and  output  data. 
The  goal  of  the  controller  is  to  continuously  estimate  the  parameters  creating  a  vector 

0c(t)  to  drive  the  error  of  the  output  vector  z(t)  between  the  reference  model  and  the 
actual  plant  to  zero. 

A  number  of  methods  for  system  modeling  and  parameter  identification,  in  terms 
of  effectiveness  and  complexity,  are  investigated  in  this  thesis.  To  improve  the  system’s 
overall  performance  and  stability  under  various  operating  conditions  and  faults,  the 
results  of  the  parameter  identification  are  applied  to  the  design  of  an  adaptive  controller 
for  a  power  grid  connected  WECS.  A  basic  view  of  the  overall  system  is  illustrated  in  the 
following  figure. 
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Since  there  are  two  independent  control  channels,  one  for  the  active  and  one  for 
the  reactive  power,  and  they  both  have  different  requirements,  the  effectiveness  of  using 
either  two  independent  adaptive  controllers  or  a  combination  of  one  adaptive  and  one 
classical  PI  controller  is  considered.  One  can  reason  that  the  reactive  power  has  simple 
control  requirements  since  it  has  to  be  kept  to  a  constant  value.  The  active  power,  on  the 
other  hand,  has  to  follow  a  profile  for  maximum  torque,  which  changes  with  wind 
velocity  and  the  aerodynamics  of  the  turbine. 

To  test  and  assess  the  effectiveness  of  the  proposed  techniques  and  their  ability  to 
adapt  to  changes  in  dynamics,  a  number  of  different  cases  have  been  simulated,  such  as 
changes  in  internal  resistances  and  wind  gusts. 

It  was  found  that  DMRAC  can  be  effectively  applied  in  a  WECS,  which  can 
significantly  improve  the  perfonnance  of  the  overall  system.  To  that  end,  a  physical 
implementation  to  validate  the  results  of  this  thesis  is  recommended. 
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I.  INTRODUCTION 


A.  BACKGROUND 

The  possibility  of  a  shortage  in  fossil  fuels  in  the  near  future  in  conjunction  with 
the  continuously  increasing  demand  for  energy  [1],  [2]  and  environmental  concerns  about 
high  emissions  has  contributed  to  a  high  interest  in  renewable  energy  (RE)  worldwide. 
Among  the  various  sources  of  RE,  one  of  the  fastest  growing  and  most  promising  is  wind 
power.  !New  wind  power  capacity  added  globally  during  2010  reached  39  GW,  which  is 
more  than  any  other  renewable  technology  and  over  three  times  higher  than  the  1 1.5  GW 
of  wind  power  added  worldwide  in  2005.  As  a  result,  existing  capacity  increased  more 
than  24%  relative  to  2009,  with  total  global  capacity  close  to  198  GW  by  the  end  of  2010 
[3].  The  growth  of  the  global  wind  power  installations  over  the  past  15  years  is 
demonstrated  in  Figure  1. 
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Figure  1.  Global  wind  power  installations,  1996  to  2010.  From  [1]. 
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In  the  U.S.,  the  wind  industry  started  to  grow  in  the  1980s;  by  1986,  the  installed 
capacity  in  California  was  close  to  90%  of  the  global  installed  wind  power.  It  later  hit  a 
plateau  during  an  electricity  restructuring  period  in  the  1990s  and  regained  momentum  in 
the  past  decade  [4]. 

The  spread  of  wind  energy  (WE)  installations  indicates  the  need  to  develop  wind 
turbines  (WT)  with  higher  efficiency.  This  involves  technology  advances  that  reduce  the 
initial  cost  of  the  WT  or  allow  the  capture  of  more  energy  during  operation,  which  results 
in  a  cost-of-energy  (COE)  decrease  [5].  One  approach  that  can  lead  to  more  effective  WT 
is  the  design  of  more  sophisticated  controllers.  The  problem  of  designing  an  adaptive 
controller  for  a  variable  speed  horizontal  axis  WT  (HAWT)  using  a  doubly-fed  induction 
generator  (DFIG)  is  addressed  in  this  thesis. 

B.  THESIS  OBJECTIVES 

The  objective  of  this  thesis  is  to  apply  various  methods  of  system  identification  to 
a  DFIG  and  investigate  the  design  of  an  adaptive  controller  for  a  power  grid  connected 
wind  energy  conversion  system  (WECS)  to  improve  the  system’s  overall  perfonnance 
and  stability  under  various  operating  conditions  and  faults.  To  achieve  those  goals, 
attempts  are  made  to  control  both  mechanical  torque  and  reactive  electric  power 
independently  using  the  rotor  voltage  components  as  a  control  input. 

C.  RELATED  WORK 

Various  researchers  over  the  past  several  years  have  addressed  the  problem  of 
efficient  control  of  WECSs.  The  complexity  of  the  overall  system  and  the  random 
conditions  under  which  WTs  operate  have  led  scientists  to  investigate  different  control 
approaches  to  improve  performance  and  stability  under  various  fault  conditions.  A  novel 
control  method  for  smoothing  the  stator  active  or  reactive  power  ripple  components  under 
unbalanced  grid  voltage  is  proposed  in  [6].  In  [7],  the  design  of  a  multi-input  multi-output 
(MIMO)  sliding  mode  (SM)  nonlinear  control  strategy  for  a  WECS  is  presented.  The 
resulting  robust  control  law  guarantees  finite  time  convergence,  whereas  smoothing 
discontinuities  at  the  basis  of  sliding  mode  control  achieves  chattering  reduction.  The 
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paper  concludes  that  this  type  of  control  is  a  particularly  suitable  option  to  deal  with 
electronically  controlled  variable  speed  operating  WECS.  Furthermore,  SM  control  has 
proven  to  be  robust  with  respect  to  system  parameter  variations  and  external  disturbances 
as  well  as  poorly  known  operating  environments  [8],  [9].  In  [10],  direct  model  reference 
adaptive  control  (DMRAC)  was  used  to  design  an  adaptive  collective  pitch  controller  for 
a  HAWT.  The  objective  of  the  adaptive  pitch  controller  was  to  regulate  generator  speed 
for  wind  speeds  higher  than  the  rated  value  for  the  generator  and  to  reject  step 
disturbances.  This  objective  was  accomplished  by  collectively  pitching  the  turbine 
blades.  The  results  were  compared  in  simulations  with  a  classical  proportional  integrator 
(PI)  collective  pitch  controller.  In  the  simulations,  the  adaptive  pitch  controller  showed 
better  speed  regulation  for  high  wind  speeds. 

A  different  approach  of  optimal  adaptive  control  is  presented  in  [11]  for  a 
permanent  magnet  synchronous  generator  driven  by  a  pulse  width  modulation  (PWM) 
inverter.  The  controller  focused  on  optimization  of  the  system  by  minimizing  Joule  losses 
in  the  inverter  and  generator.  To  achieve  this,  the  paper  proposed  the  operation  of  the 
system  at  an  operating  point  with  minimum  phase  current,  considering  fixed  torque. 

Other,  more  classical  approaches  include  the  independent  control  of  active  and 
reactive  power  using  back-to-back  PWM  voltage  source  converters  in  the  rotor  circuit.  A 
vector  control  scheme  on  the  supply  side  PWM  ensures  independent  control  of  active  and 
reactive  power  drawn  from  the  supply,  while  vector  control  of  the  rotor  side  PWM  results 
in  a  wide  speed  -  range  operation.  The  control  is  achieved  using  PI  controllers  to  enable 
maximum  speed  tracking  for  maximum  power  extraction  from  the  wind.  This  approach, 
presented  in  [12],  presents  two  different  controls,  a  current-mode  control  and  a  speed¬ 
mode  control.  A  similar  approach  was  presented  in  [13],  where  a  Stator-Flux  Oriented 
vector  control  approach  is  deployed  for  both  stator  and  rotor  side  converters  to  provide 
independent  control  of  active  and  reactive  power  using  a  decoupled  design  based  on 
internal  model  control. 
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D.  APPROACH 


In  the  first  part  of  this  research,  the  analytical  design  of  an  online  identification 
algorithm  is  presented,  with  the  intent  of  constructing  a  Simulink  model  of  the  physical 
system.  In  the  second  part,  an  adaptive  control  scheme  to  compensate  for  dynamic 
variations  of  the  generator  is  proposed  and  its  effectiveness  is  evaluated  by  computer 
simulations.  In  the  simulations,  various  fault  conditions  are  implemented,  and  the  results 
are  discussed  subsequently. 

E.  THESIS  OUTLINE 

This  thesis  is  organized  as  follows:  the  concept  of  the  basic  theory  of  WTs  and 
DFIGs  is  covered  in  Chapter  II,  and  the  model  created  for  the  system  in  the  Simulink 
environment  is  described.  The  theory  of  online  system  identification  including  modeling 
and  simulation  results  for  different  estimation  algorithms,  is  covered  in  Chapter  III.  A 
basic  introduction  to  adaptive  control  theory  and  the  design  and  simulation  of  two 
different  adaptive  controllers,  along  with  simulation  results  under  different  operational 
conditions  of  the  system,  are  presented  in  Chapter  IV.  A  comparison  with  simulation 
results  from  other  control  approaches  is  also  presented.  Finally,  the  conclusions  based  on 
the  results  from  the  previous  chapters,  along  with  some  recommendations  for  future 
research  in  the  area  of  adaptive  control  for  WTs,  are  presented  in  Chapter  V.  The 
Appendix  includes  Matlab  code  for  initialization  of  the  simulation  of  the  controllers  as 
well  as  the  basic  model  blocks  in  the  Simulink  environment.il.  WIND  TURBINES  - 
DFIG 

A.  WIND  TURBINE  THEORY 

Wind  turbines  are  designed  to  capture  a  portion  of  the  wind's  kinetic  energy  and 
convert  it  to  usable  energy  [14]  in  the  form  of  electricity.  This  section  presents  a  brief 
explanation  of  the  basic  operating  principles  of  WTs. 
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1.  Energy  in  the  Wind 


Wind  speed  is  a  random  process  that  depends  on  several  factors.  The  main 
parameter  that  describes  the  ability  to  extract  power  is  the  mean  wind  speed  for  the 
specific  area  of  interest.  The  probability  distribution  function  (PDF)  of  the  mean  wind 
speed  Vm  over  an  area  is  modeled  as  a  Weibull  distribution  [14]  of  the  form: 
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where  the  parameters  k  and  c  are  called  the  shape  and  scale  factor,  respectively,  and  are 
both  positive  real  numbers. 
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Figure  2.  A  volume  of  air  flowing  with  speed  V. 


If  one  takes  a  volume  of  air  AAx  flowing  with  speed  V  as  shown  in  Figure  2,  the 
kinetic  energy  stored  in  the  wind  is 


E 


k 


(0.0) 


where  m  is  the  mass  of  the  air. 

If  we  replace  the  mass  in  terms  of  the  air  density  p  and  the  volume  as  m  -  pAAx , 
Equation  (2.2)  can  be  written  as 
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Ek=^p(AAx)V2 , 

and  from  Equation  (2.3)  one  can  detennine  the  wind  energy  per  unit  volume  AAx  =  1  as 

E,^pV\  (0.0) 

Since  the  power  is  defined  as  the  derivative  of  the  energy  with  respect  to  time, 
one  can  derive  an  expression  for  the  power  of  the  wind: 


AE 

At  At 


(0-0) 


2.  Wind  Turbine  Aerodynamics  -  Disc  Actuator  Model 

Turbine  aerodynamics  describes  the  forces  developed  on  a  wind  turbine  by  an 
airflow  passing  through  it.  The  approach  used  to  derive  a  model  for  the  wind  turbine  is 
actuator  disc  theory  [14].  In  this  theory,  the  WT  is  regarded  as  an  actuator  disc,  a  device 
that  is  designed  to  extract  part  of  the  kinetic  energy  of  the  wind.  This  theory  also  gives  a 
value  for  the  maximum  efficiency  limit.  If  one  assumes  the  disc  is  immersed  in  an  airflow 
as  shown  in  Figure  3,  then  the  upstream  wind  speed  V  has  to  be  greater  than  the 
downstream  wind  speed  V-x. 
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A.oo 


Furthermore,  for  the  stream  tube  just  enclosing  the  disc,  the  upstream  cross 
sectional  area  is  smaller  than  the  disc  area  Ad,  which  in  turn  is  smaller  than  the 
downstream  cross-sectional  area  A. This  happens  because,  by  definition,  due  to  the  law 
of  mass  conservation,  the  mass  flow  must  remain  constant  within  the  stream  tube  [14]. 
The  above  statement  implies: 


AM 

- =  constant , 

At 


(0.0) 


where  AM  is  the  air  mass  going  through  a  surface  during  time  At.  Now  one  can  derive  the 
following  expression: 


=  pA  —  =  pAV  =  constant.  (0.0) 

At  At 


Therefore,  if  one  takes  ideal  points  away  from  the  turbine  (-  and  +  infinity)  and 
on  the  turbine,  one  obtains: 


pAy  =  PADvD  =  pA_y_m 


(0.0) 


7 


Since  the  air  that  passes  through  the  disc  undergoes  a  speed  reduction  equal  to 
V -V_x,  this  implies  a  reduction  in  its  kinetic  energy.  This  results  in  a  force  FD  developed 

by  the  actuator  disc  on  the  incident  airflow,  which  is  given  by  Newton's  second  law  of 
motion  as  the  product  of  mass  M  and  acceleration  a: 

Fd  =  pA0 Ax^-^  =  (V  -  V_JpADVD .  (0.0) 

At 

The  force  FD  originated  by  the  pressure  drop  introduced  by  the  actuator  disc  can 
also  be  expressed  as: 

Fd=(Pd-Pd)ad »  (0.0) 

where,  pD  and  pD  are  the  values  of  the  air  pressure  immediately  before  and  after  the  disc, 
as  shown  in  Figure  4. 


Figure  4.  Pressure  around  actuator  disk.  After  [15]. 

Since  Vd  is  smaller  than  V,  one  can  write  the  speed  at  the  disc  as 

VD=(1  -a)V,  (0.0) 

where  one  defines  a,  with  0  <  a  <  1 ,  which  is  referred  to  as  the  “axial  flow  interference 
factor.”  Substituting  Equation  (2.1 1)  in  Equation  (2.9),  one  gets  the  force  FD  as: 

FD=(V-V_JpADV(l-a).  (0.0) 
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Equating  Equations  (2.10)  and  (2.12),  one  obtains: 

(V-VjpV(l-a)  =  (p+D-pD).  (0.0) 

For  an  airflow  that  is  stationary,  incompressible,  and  frictionless  and  with  no 
external  forces  applied  either  upstream  or  downstream,  the  Bernoulli  equation  states  that 
the  total  energy  of  the  airflow  remains  constant.  Given  that  the  above  conditions  are 
satisfied,  one  can  apply  the  Bernoulli  equation  from  far  upstream  to  just  in  front  of  the 
rotor  and  from  just  behind  the  rotor  to  far  downstream  [16]: 

(p'o-Po)  =  \p{V2-Vl).  (0.0) 

Substituting  Equation  (2.14)  in  Equation  (2.13),  one  finally  gets 

l/00=(l-2  a)V,  (0.0) 

and,  since  V  x  has  to  be  non-negative,  Equation  (2.15)  gives  us  the  upper  bound  of  a, 
which  is  0.5. 

Substituting  Equation  (2.15)  in  Equation  (2.12),  one  gets  an  expression  for  the 
force  that  is  applied  from  the  airflow  on  the  disc: 

FD=2pADV2a(l-a).  (0.0) 

The  power  that  the  actuator  disc  extracts  from  the  wind  is  given  by 

PD=FDVD=2pADV3a(l-a)2.  (0.0) 

This  gives  an  idea  of  the  maximum  power  available  to  the  turbine  from  the  flow 
of  the  wind. 

A  conventional  way  of  characterizing  the  ability  of  a  wind  turbine  to  capture  wind 
energy  is  the  power  coefficient  Cp  [14],  which  is  defined  as  the  ratio  of  extracted  power 
divided  by  the  total  wind  power: 
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Substituting  Equations  (2.2)  and  (2.11)  in  (2.12),  one  gets  the  following 
expression  for  Cp: 


C 


P 


2pADV3a(l-a)2 
0.5  pADV3 


=  4a(l-a)2 . 


(0.0) 


The  power  coefficient  Cp  is  a  function  of  the  axial  flow  interference  factor  a.  The 
maximum  theoretical  value  of  the  power  coefficient  can  be  found  from  Equation  (2.19) 
by  taking  its  first  and  second  derivatives.  The  result,  which  is  called  the  Betz  limit,  is 
Cp  max  =  0-593  for  a  equal  to  1/3  [17].  This  limit  is  valid  for  all  kinds  of  wind  turbines,  but 

the  actual  efficiency  is  usually  lower,  and  for  most  commercial  wind  turbines,  a  good 
value  is  between  0.40-0.50. 

In  practice  three  effects  that  lead  to  a  decrease  in  the  maximum  achievable 
power  coefficient  [18]: 

•  Rotation  of  the  wake  behind  the  rotor. 

•  Finite  number  of  blades  and  associated  tip  losses. 

•  Non  zero  aerodynamic  drag. 

3.  Torque  and  Power  on  the  Rotor 

The  relation  between  the  energy  of  the  wind  flowing  through  the  windmill  and  the 
mechanical  power  available  at  the  rotor  is  discussed  in  this  section.  The  parameter  [3  is 
the  pitch  angle  of  the  turbine  blades,  and  X  the  tip  speed  to  wind  speed  ratio  : 

A  =  ^~.  (0.0) 

V 


A  graphical  representation  of  the  parameters  involved  in  Equation  (2.20)  is  shown 
in  Figure  5. 
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Figure  5.  Turbine  disc. 


Substituting  the  expression  for  the  power  coefficient  Cp  from  Equation  (2.19)  into 
Equation  (2.17),  one  gets  the  following  expression  for  the  aerodynamic  wind  power: 

PD  =2pAyC^^-  =  ±pADCp(X,P)V\  (0.0) 

Substituting  the  expression  for  the  disc  area  AD  -  ttR2  into  Equation  (2.21),  one 


gets: 

PD=^p7iR2Cp(A,j3)V3.  (0.0) 

From  Equation  (2.22),  one  can  easily  derive  an  expression  for  the  aerodynamic 
torque  produced  by  a  turbine  of  radius  R  rotating  at  an  angular  rate  Q,  with  wind  speed  V: 


T  = 


Pr  _1  pxR2Cp(A,j3)V3 


(0.0) 


2  Qr 

In  Equation  (2.23),  the  substitution  of  Qr  by  its  value  extracted  from  Equation 
(2.20)  provides  the  following  expression: 

Tr=^P7rR3CQ(A,j3)V2 

where  the  Cq  is  named  the  torque  coefficient  and  is  related  to  Cp  by 


(0.0) 
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In  Equations  (2.21)  and  (2.22),  it  can  be  seen  that  Cq  and  Cp  are  written  as 
functions  of  X  and  fi.  However,  it  ia  assumed  in  this  thesis  that  the  pitch  angle  /?=(),  and 
the  above  coefficients  are  used  as  functions  of  the  ratio  X  only. 


A  parametric  approximation  of  the  torque  coefficient  Cq  is  given  by  the  following 
expression  in  terms  of  X  [7] 


ceW  =  j 


*-l] 

2  J 


-c!  X 


(0.0) 


where  a,  b,  c  are  dimensionless  parameters  that  characterize  the  rotor  of  the  wind  turbine. 
If  one  substitutes  Equation  (2.26)  in  Equation  (2.25),  then 


Cp(X)  =  a 


2  1 


—c!X 


(0.0) 


The  plot  of  a  typical  power  coefficient  curve  versus  X  is  shown  in  Figure  6. 


Figure  6.  Typical  Cp  versus  speed  ratio  X  curve. 


The  speed  ratio  X  is  a  very  important  parameter  in  the  dynamic  model  of  the  WT. 
Since  one  wants  to  operate  as  closely  as  possible  to  maximum  efficiency  to  have 
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maximum  power  capture  and  given  that  the  wind  speed  cannot  be  controlled  and  the 
radius  of  the  blades  is  fixed,  one  needs  to  vary  the  rotational  speed  of  the  wind  turbine  to 
keep  /l  close  to  its  optimal  value  Xopt.  This  is  one  of  the  advantages  of  variable  speed  WTs 
over  fixed  speed  WTs,  since  they  operate  close  to  optimal  efficiency  only  for  a  small 
range  of  wind  speed  values. 

For  a  typical  wind  speed  profile  with  a  Weibull  distribution  having  a  shape 
parameter  k  =  2  and  scale  parameter  c  =  8.5,  the  variable-speed  turbine  captures  2.3% 
more  energy  per  year  than  the  constant-speed  turbine.  This  is  considered  a  significant 
difference  in  the  wind  industry  [18]. 

4.  Regions  of  Control 

There  are  three  different  operational  regions  for  the  variable  speed  WTs,  which 
are  shown  in  Figure  7. 
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In  region  1,  the  section  of  the  solid  curve,  usually  located  to  the  left  of  4  to  5  m/s 
wind  speed,  the  power  available  in  the  wind  is  low  compared  to  the  losses  in  the  turbine 
system;  therefore,  the  turbine  is  usually  not  operating  in  this  region.  The  use  of  modem 
control  strategies  is  not  usually  critical  in  region  1  [5]. 

Region  2,  illustrated  by  the  cubic  section  of  the  red  curve,  is  generally  for  winds 
speeds  between  4  to  5  m/s  up  to  the  WT’s  rated  wind  speed,  usually  between  12  to  15 
m/s.  In  this  area,  the  WT  is  in  operational  mode,  and  it  is  desirable  for  it  to  capture  as 
much  power  from  the  wind  as  possible.  The  goal  in  this  region  is  to  approach  the  Betz 
Limit  curve  as  closely  as  possible. 

Region  3  is  for  operation  above  rated  wind  speed  of  the  WT,  i.e.,  the  wind  speed 
above  which  maximum  peak  power  is  produced.  The  goal  of  control  strategies  applied  in 
this  region  is  to  limit  the  wind  power  captured  so  as  not  to  exceed  the  designed  electrical 
and  mechanical  load  limits  of  the  system. 

All  research  in  this  thesis  centers  on  Region  2. 

B.  DOUBLY-FED  INDUCTION  GENERATOR  (DFIG) 

Induction  machines  are  used  in  a  wide  variety  of  applications  and  are,  without 
doubt,  the  workhorse  of  the  electric  power  industry  [19]. 

An  induction  motor  is  an  electric  motor  that  is  driven  by  alternating  current  (AC). 
It  consists  of  two  basic  parts:  the  stator,  which  is  on  the  outside  and  has  windings 
supplied  with  AC  to  produce  a  rotating  magnetic  field,  and  the  rotor,  which  is  attached  to 
the  output  shaft  on  the  inside.  When  balanced  multi-phase  currents  flow  through  the 
stator  windings,  a  rotating  air-gap  magnetomotive  force  (MMF)  is  produced,  and,  if  the 
rotor  rotates  at  different  speed  than  the  speed  of  the  rotating  MMF,  balanced  multi-phase 
currents  are  induced  in  the  rotor  windings.  There  are  two  different  types  of  rotor 
windings:  wound  and  squirrel  cage.  In  Figure  8  the  cross  section  of  a  typical  three-phase 
induction  machine  with  a  squirrel  cage  rotor  is  illustrated. 
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Figure  8.  Induction  machine  cross  section.  From  [20]. 


In  single-fed  induction  generators,  the  rotor  windings  are  short-circuited,  while 
DFIGs  take  full  advantage  of  both  stator  and  rotor  excitations.  The  rotor  currents  and 
voltages  of  the  DFIG  are  at  the  basis  of  the  control  techniques  developed  in  this  thesis. 

1.  Reference  Frames  -  Dynamic  Models 

Three  phase  signals  (currents  or  voltages)  can  be  represented  either  by  the 
respective  time  domain  components  (called  “abc”)  or,  more  conveniently,  by  the 
expansion  with  respect  to  a  rotating  reference  frame  (called  “qdO”). 

The  equations  converting  time  domain  "abc"  voltages  into  "qdO"  components  and 
vice  versa  are  very  standard  and  can  be  found  in  a  number  of  sources  [12],  [19]. 

In  the  "qdO"  reference  frame  for  a  symmetrical  DFIG,  we  have  the  following 
equations  for  the  flux  linkages  from  [19]: 
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PXqd Os  =-rJqdOs  -j<*>AdqOs+VqdOs 
PKdor  =  +  j(& ~  KrfOr  +  V'qdOr 

where  v,  i,  X  refer  to  voltages,  current,  and  flux  linkages  in  the  qdO  reference  frame  of  the 
stator  (subscript  5)  and  the  rotor  (subscript  r).  The  parameter  co  is  the  electrical  angular 
rate  of  the  reference  frame,  determined  by  the  line  frequency  (say  60  Hz  in  the  USA);  cor 
is  the  electrical  frequency  of  the  rotor  circuit;  rs  ,rr  are  the  stator  and  rotor  resistances; 
and  p  stands  for  the  differential  operator  d/dt. 

Any  qdO  variable,  say /,  can  be  expressed  in  a  complex  form  as: 

fqdO=fq~jfd+0  ■  (0-0) 


Therefore,  the  flux  linkages  are  defined  as: 

(4,os)T=[4  -4s  0] 

(4U)r=K  -4,  o]' 


(0.0) 


where  the  flux  components  are  related  to  the  currents  by  standard  equations  [19] 


4s  -LlJqs  +LM('qs+'qr ) 

4s  =  LJds  +  Lm  [ids  +  1 dr  ) 
4s  =  i-Jos 

4r  =4V  +LM('qs  +'qr ) 
4r  =  i-lddr  +  [' ds  +  'dr  ) 
4r  =i-lr'or 


(0.0) 


Expanding  Equations  (2.28)  and  combining  them  with  Equations  (2.31),  one  gets 
the  following  set  of  equations  to  characterize  the  DFIG  [19]: 


16 


rs+^*ss 

CO  .. 


A* 


r+P  X„ 


G>b 

0 


P_x 

- 


Ay, 


-Ay, 


r+Ay 


Ax 


Ax 

—  am 
P  V 


f  CO  -  6)r  ) 

X 

r+  AX  1 

'  co  -  cor  ) 

l  J 

0 

/  \ 

V  J 

r+Ax’ 

0 


0 

0 

r>  Ax„ 


(0.0) 


where 


Xss  =  X,s+XM 


xrr  =  xlr  +  xM 

X  Is  —  ^b^ls 


(0.0) 


X/r  (0b^!r 


X  M  ^b^M 


express  all  the  impedances.  In  the  above  equations,  cob  is  the  base  or  rated  frequency  (say 
the  standard  60  Hz);  cor  is  the  electrical  rotational  frequency  of  the  rotor;  and  coe  is  the 
electrical  rotational  frequency  of  the  stator.  The  variables  Lh  and  LM  are  the  leakage  and 
magnetizing  inductances  of  the  stator  windings,  and  L  h  is  the  leakage  inductance  for  the 
rotor  windings.  Also,  as  mentioned,  rs  and  r',  are  the  stator  and  the  rotor  resistances. 


The  electromagnetic  torque  of  the  generator  is  given  by  the  following  expressions 
[19],  which  are  all  equivalent: 


T  = 


( 31 

f  PJ 

UJ 

UJ 

f'/  -  lyi;'  ^  5 


(0.0) 


T  = 


T  = 


(  3^| 

( u 

UJ 

UJ 

f  3^| 

UJ 

UJ 

Xdriqr) 


^-M^XdJqs  XqJdi) 


(0.0) 


In  the  above  expressions,  P  is  the  number  of  magnetic  poles  of  the  generator,  and 
Te  is  positive  for  motor  action  and  negative  for  generator  action. 
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The  total  real  and  reactive  power  of  a  three-phase  balanced  DFIG  in  steady  state 
conditions  is  given  by  [19]: 


Pe=\wqJqs  +  vdjds) 

a=fov*-W 


(0.0) 


2.  DFIG  Model 


Equations  (2.28)  clearly  define  a  nonlinear  dynamic  system  where  the  fluxes  X  are 
the  states  and  the  voltages  are  the  control  signals.  Equivalently,  from  [19],  one  can 
redefine  the  states  in  terms  of  the  qdO  current  components  of  the  stator  and  the  rotor. 

The  whole  system  is  a  nonlinear  dynamic  system  and  was  implemented  by 
Professor  A.  Julian  in  Simulink  [21],  as  shown  in  Figure  9.  The  time  domain  voltage 
components  va,  v&  vc  of  the  stator  and  the  rotor  are  explicitly  shown,  together  with  the 
state  variables  ia,  ib,  4  of  stator  and  rotor,  in  the  time  domain.  In  addition,  other  physical 
quantities,  like  the  torque  on  the  load  TLi  are  shown,  along  with  all  relevant  power 
components  Pe  (active)  and  Qe  (reactive).  The  model  of  the  DFIG  is  described  in  more 
detail  in  [21]. 


Figure  9.  General  view  of  the  DFIG  model. 
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This  model  is  used  extensively  in  the  next  chapters,  both  as  a  benchmark  for 
simulating  the  dynamics  of  the  DFIG  and  as  the  basis  of  the  analytical  design  of  the 
proposed  algorithms. 

C.  CHAPTER  SUMMARY 

The  basic  theory  behind  WTs  and  DFIGs  was  presented  in  this  chapter,  and  a 
description  of  the  simulation  model  for  a  DFIG  was  provided.  The  quality  of  the  model  is 
of  great  importance,  because  the  effectiveness  of  the  controller  depends  on  it.  In  the  next 
chapter,  three  different  algorithms  of  online  parameter  identification  will  be  presented. 
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III.  ONLINE  PARAMETER  ESTIMATION 


A.  INTRODUCTION 

The  mathematical  model  described  in  the  Chapter  II,  with  its  Simulink 
implementation,  is  representative  of  the  DFIG  physical  system.  In  general,  the  various 
parameters  defined  (resistances  and  inductances)  have  to  be  estimated  based  on  voltage 
and  the  current  measurements  of  the  physical  system  itself.  The  mathematical  model  can 
be  used  for  control  design  and  validation  by  computer  simulation. 

The  purpose  of  this  chapter  is  to  describe  the  methods  used  to  achieve  on-line 
estimation  of  the  DFIG  parameters.  The  essential  idea  behind  on-line  estimation  is  the 
comparison  of  the  observed  system  response  y(t)  with  the  output  of  a  parameterized 
model  whose  structure  is  the  same  as  that  of  the  plant  model  [22].  In  this  case,  the  output 
vector  y(t)  consists  of  all  the  measured  currents  in  the  qdO  reference  frame  under  different 
conditions  of  excitation  through  the  input  voltages. 

The  parameter  vector  6{t)  contains  estimates  of  all  the  parameters  in  the  state 
space  model  in  Chapter  II,  and  it  is  adjusted  continuously  so  that 
y[6,t)  approaches  y(0  as  t  increases.  Under  certain  input  conditions,  y(0,t)  being  close 

to  y(t)  implies  that  ()(t)  is  close  to  the  unknown  parameter  vector  9(t)  of  the  plant  model. 
The  online  estimation  procedure,  therefore,  involves  three  steps  [22]: 

•  Create  an  appropriate  parameterization  of  the  plant  model; 

•  Design  an  algorithm  for  generating  and  updating  the  parameter  vector 
estimate;  and 

•  Properly  excite  the  plant  such  that  the  parameter  estimates  approach  the 
actual  expected  values  as  t  — »  oo  . 

In  general,  the  plant  model  of  a  system  to  be  studied  for  identification  should  be 
expressed  in  the  following  fonn: 
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(0.0) 


y_(t)  =  0T(tm 


which  can  be  explicitly  written  as: 


kl(t) 

y2(t) 

<Pu(t)  <Pu(t)  . 
(p21[t) 

■  ■ 

m 

yn(t)_ 

where  the  nx  1  vector  v(t)  consists  of  the  observed  variables;  the  mx  1  vector  9(t)  consists 
of  the  unknown  parameters  to  be  determined;  and  the  nxm  matrix  (PT(t)  consists  of  known 
functions  that  describe  the  plant  model,  which  are  called  regression  variables[23],  In  this 
thesis,  as  one  can  see  later  on,  n=6  and  m= 5. 

It  is  shown  in  [22],  [23]  that  under  conditions  of  persistency  of  excitation,  the 
estimated  parameters  converge  to  the  true  parameters.  Loosely  speaking,  this  is  satisfied 
when  the  signal  is  sufficiently  excited  so  that  the  input-output  behavior  can  be  modeled 
by  a  unique  parameter  vector. 

B.  SYSTEM  EQUATIONS 

Based  on  what  was  discussed  in  the  previous  section,  one  needs  to  rewrite  the 
dynamic  equations  in  terms  of  parameters  and  regression  vectors.  In  particular,  the  five 
quantities  that  form  the  parameter  vector  9  to  be  estimated  are  the  following: 

0  =  \rs  r'r  Lls  C  Lm]T  ■  (°-°) 

The  vector  y(t)  of  observations  is  defined  by  the  stator  and  rotor  voltages: 

y(t)  =  v(t)  =  [vqs  vds  v0s  vqr  vdr  v0r]T .  (0.0) 

Substituting  Equations  (3.3)  and  (3.4)  in  Equation  (2.32),  one  obtains: 
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V(t)  = 


'Os 

0 

0 

0 


0 

0 

0 


Piqs+ojids 
P'ds  -  0J'qs 

pi  os 
0 
0 
0 


0 

0 

0 

Piqr  +(CO~CDr)idr 

pidr-(co-o)r)iqr 

PI  Or 


p(iqs+iqr)  +  o)(ids+idr) 

P(ids+idr)-^(iqs+iqr) 

0 

P(iqs+iqr)  +  (v-Mr)(ids+idr) 
P(ids+idr)-(v-Mr)(iqs+iqr) 
0 


0(t),  (0.0) 


where  the  matrix  on  the  right  of  Equation  (3.5)  becomes &T  (t)  in  Equation(3.2): 


V(t)=ojmt), 


(0.0) 


and  co  is  the  angular  velocity  of  the  arbitrary  reference  frame  that  was  chosen. 

C.  ESTIMATORS 

In  order  to  estimate  the  machine  parameters  vector  9(t),  three  different 
identification  approaches  were  used: 

•  least  mean  squares  (LMS), 

•  recursive  least  Squares  (RLS),  and 

•  recursive  least  squares  with  exponential  forgetting  (RLS-EF). 

The  next  section  describes  the  three  different  estimators. 

1.  Least  Mean  Squares  (LMS) 

The  first  estimation  algorithm  is  based  on  the  Least  Mean  Squares  (LMS)  method. 
This  was  initially  formulated  at  the  end  of  the  18th  century  by  the  German  mathematician 
and  scientist  K.  F.  Gauss  to  calculate  the  orbits  of  planets  and  asteroids.  From  [23],  using 
the  above  method,  the  parameters  of  the  model  should  be  chosen  in  such  a  way  that  the 
sum  of  squares  of  the  differences  between  the  observed  and  computed  values,  multiplied 
by  numbers  that  measure  the  degree  of  precision,  is  a  minimum. 
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a. 


Equations  Setup 


Based  on  the  model  described  in  Equation  (3.1),  one  can  define  the 
prediction  or  estimation  error  eft)  [23],  as  the  difference  between  y(t)  and  y(t ) : 

e(t)  =  y(t)-y(t)  =  y(t)-OT(t)0(t).  (0.0) 

In  addition,  one  can  define  the  parameter  error  Oft)  as  follows: 

(0.0) 

Thus,  from  Equations  (3.1),  (3.7)  and  (3.8),  one  can  derive  the  following 
expression  for  eft): 


e(t)=^mr) 


(0.0) 


From  [24],  the  general  form  of  an  identification  algorithm  is  defined  by  a 
differential  equation,  called  the  update  law,  which  is  of  the  form: 

S(t)  =  F{y(t),e(t),mMt)}.  (0.0) 


By  the  standard  steepest  descent  algorithm,  one  can  obtain  the  following 
expression  for  the  derivative  of  Oft)  with  respect  to  time: 

o(t)=-<u®mt),  (o.o) 


where  u  is  an  arbitrarily  chosen  strictly  positive  parameter  called  the  adaptation  gain  and 
allows  for  a  variable  rate  of  adaptation  of  the  parameters.  The  greater  the  value  of  ju,  the 
faster  the  convergence,  but  the  response  is  more  sensitive  to  noise  and  might  induce 
oscillations  in  the  system. 


function: 


To  analyze  the  behavior  of  this  system,  one  can  define  a  Lyapunov 


v(0(t))  =  -  Oft) 


(0.0) 
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According  to  [25],  for  0[t)  to  converge,  the  Lyapunov  function  described 
in  Equation  (3.12)  has  to  be  positive  definite,  and  its  total  derivative  with  respect  to  time 
must  be  negative  semi-definite.  By  definition,  it  is  positive  definite,  since  in  this  case  the 
parameters  are  positive  real  numbers  and  this  case  takes  the  norm  of  the  parameter  vector. 
The  derivative  with  respect  to  time  of  V(0(t))  is  given  by: 

mt))=^^l=dT(t)d(t).  (0.0) 

dt 

Substituting  the  expression  of  6[t)  from  Equation  (3.11)  into  Equation 
(3.13),  one  can  obtain: 

V(0,t)  =  -judT(tm)e(t).  (0.0) 

Since  the  prediction  error  e(t)is  related  to  the  parameter  error  by 
Equation  (3.9),  one  can  obtain: 

V(0,t)  =  -//£T(t)®(t)Or(t)£(t)  =  -JU I  ®T(t)m  |2<  0  .  (0.0) 

The  outcome  of  Equation  (3.15)  is  that  V(0,t)  is  negative  semi-definite, 
and  as  a  result,  6[t)  decreases  monotonically  with  time.  Under  conditions  of  persistency 
of  excitation,  this  implies  that  the  parameter  error  actually  converges  to  zero  as: 

0[t)~kt)->  0.  (0.0) 

Since  the  system  parameter  vector  0{t)  is  assumed  constant  or  slowly 
changing  with  time,  one  can  assume  that  0{t ) «  0 .  Therefore,  by  differentiating  Equation 
(3.8)  and  substituting  Equation  (3.11),  one  gets  the  following  expression  for  the 
derivative  of  the  parameter  estimates: 

d(t)  =  -0(t)  =  v®(t)e(t).  (0.0) 

Integrating  Equation  (3.17),  one  gets  the  expression  needed  for  the 

estimates: 
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(0.0) 


m  =  m  +  J  jU®(T)[y(T)  -  ®T(r)£(r)]dr , 

0 

where  the  initial  condition  can  be  arbitrarily  chosen  based  on  any  a  priori  knowledge  of 
the  plant. 

b.  Implementation 

The  simulation  of  this  estimation  method  required  the  creation  of  a 
Simulink  subsystem  which  was  added  to  the  model  described  in  Chapter  II  for  the  DFIG. 
A  general  view  of  the  estimator  is  shown  in  Figure  10. 
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Figure  10.  General  view  of  the  LMS  estimation  model. 


The  Simulink  blocks  that  provide  the  vectors  0(t)  and  V(t)  are  common 
for  all  the  different  methods  and  were  implemented  using  Equations  (3.6).  These  are  the 
equations  for  reference  frame  transfonnation  described  in  [19].  The  estimator  block  for 
the  LMS  is  shown  in  more  detail  in  Figure  1 1 . 
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Figure  1 1 .  Detailed  block  diagram  for  the  LMS  estimator  model. 


It  can  be  seen  in  Figure  1 1  that  instead  of  using  the  same  adaptation  gain  // 
for  all  the  parameters,  different  values  of  //  for  the  different  channels  were  employed.  The 
adaptation  matrix  in  Figure  1 1  is  a  diagonal  matrix  of  the  form: 


/u_Matrix 


fin  0  0  0  0 

0  jiu22  0  0  0 

0  0  ju33  0  0 

0  0  0  //44  0 

0  0  0  0  ju55 


The  reason  behind  this  is  to  compensate  for  different  scaling  factors  and 
bring  all  the  signals  to  a  similar  order  of  magnitude.  This  seems  to  improve  the 
convergence  rate  and  prevent  instability. 

The  values  of  the  non-zero  terms  in  the  //  matrix  were  chosen  to  create  a 
relatively  fast  convergence  without  affecting  the  stability  of  the  system.  Large  values  of 
the  parameters  jin  lead  to  fast  convergence,  but  the  estimates  are  oscillating.  Very  small 
values  lead  to  very  slow  but  stable  convergence.  To  illustrate,  the  values  chosen  for  the 
above  parameters  that  gave  a  stable  output  are  shown  in  Table  1. 
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Table  1.  Adaptation  gains  for  the  LMS  simulation. 


Pii 

Value 

Pll 

130 

P22 

820 

1133 

0.007 

P44 

0.38 

P55 

0.0135 

2.  Recursive  Least  Squares  (RLS) 

The  second  method  was  recursive  least  squares  (RLS). 


a.  Equations  Setup 


This  method  is  implemented  in  discrete  time  and  recursively  computes  the 
parameter  estimates  9[n)  by  minimizing  the  quadratic  error  [23]: 


V(9,n)  =  ^fj[y(iT5)-0T(iT5)9]2, 

2  /=! 


(0.0) 


with  Ts  as  the  sampling  interval. 

From  [26],  the  recursive  estimate  can  be  interpreted  as  a  Kalman  Filter  for 

the  process: 


9[n  + 1)  =  9(n) 
y(n)  =  <&J(n)9(n)  +  e(n) 


(0.0) 


where  y(n)  <P(n)  are  given  by  Equation  (3.4)  and  Equation  (3.6),  respectively,  by 
substituting  t  with  n. 

The  recursive  equations  for  the  RLS  can  then  be  derived  from  the  Kalman 
fdter  equations  used  in  [27]: 
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(0.0) 


6{n  + 1)  =  Q[n)  +  K(n)[y(n  + 1)  -  0>  +  l)d(n)] 

K[n)  =  P(n)0(n  +  l)[0T(n  +  l)P(n)0(n  + 1)  +  /]  1 
P(n  + 1)  =  P(n)  -  /C(n)[0T(n  +  l)P(n)0(n  + 1)  +  /]  1/Cr(n) 

where  the  matrix  ®(n)  has  to  be  full  rank,  which  means  that  it  has  to  be  nonsingular  for 
all  n  >  0 .  The  initial  condition  for  Equations  (3.21)  is: 

rn=cr2imXm  (o.o) 

where  a  is  a  positive  number,  much  higher  than  the  parameters,  and  m  is  the  dimension  of 
y(t),  which  in  this  case  is  six. 

b.  Implementation 

The  general  outlook  of  the  RLS  estimator  is  shown  in  Figure  12. 


TIME  DELAY 


Figure  12.  General  view  of  RLS  estimator  model. 


A  more  detailed  view  of  the  inside  of  the  RLS  block  is  shown  in  Figure 
13.  The  inputs  to  the  block  are  ®(n)  and  v(n),  and  the  output  is  the  parameter  vector 

estimates  + 
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Figure  13.  Detailed  view  of  RLS  estimator  model. 


3.  RLS  With  Exponential  Forgetting  (RLS-EF) 

In  the  above  estimations,  the  parameters  are  assumed  constant  or  slowly  changing. 
Now  consider  the  case  of  varying  parameters  from  either  environmental  or  operational 
conditions  or  a  change  in  the  parameters  due  to  a  fault  in  the  system.  The  case  that  is  of 
more  interest  is  when  the  parameters  are  slowly  changing  in  time.  In  this  case,  one  needs 
to  design  an  algorithm  with  a  “forgetting  factor”  which  discounts  data  exponentially 
further  back  in  time. 


a.  Equations  Setup 

According  to  [23],  the  RLS-EF  algorithm  is  obtained  by  changing  the  cost 
function  of  Equation  (3.19)  from  the  one  that  was  used  for  the  RLS  algorithm  to  the 
following: 

V(0,t)  =  ^fjAt-i[y(i)-OTmi)]2,  (0.0) 

^  i=1 
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where  X  is  a  parameter  such  that  0  <  A  <  1  and  is  called  the  forgetting  factor.  The  meaning 
of  Equation  (3.23)  is  that  one  weights  the  last  data  by  1,  with  a  factor  of  Xn  for  data  that 
are  n  time  units  old.  This  means  that  one  reduces  exponentially  the  weight  of  each  set  of 
measurements.  Thus,  the  equation  for  the  RLS  with  exponential  forgetting  becomes: 


6(t  + 1)  =  6(t)  +  K  (t)[y(t  + 1)  -  ®T(t  + 1  )0(t)] 

K(t)  =  ^®(t  +  l)[0T(t  +  + 1)  +  /]  1 

P(t  + 1)  =  —  -  K(t)[0J(t  + 1) — 0(t  + 1)  +  /]  1  Kt  (t) 


A 


A 


(0.0) 


b.  Implementation 


The  general  outlook  of  the  RLS-EF  estimator  is  similar  to  that  shown  in 

Figure  12. 


D.  SIMULATION  OF  THE  ESTIMATORS 


1.  Equipment 

The  three  estimators  described  in  the  previous  pages  were  simulated  using 
Simulink  software.  A  Matlab  code  was  created  to  initialize  the  models.  The  code  used  for 
initialization  of  the  simulation  models  in  Simulink  is  presented  in  the  Appendix.  To 
determine  how  well  the  estimators  perform,  one  has  to  choose  a  DFIG  with  known 
parameters  for  comparison.  The  DFIG  that  was  used  is  the  Lab- Volt  model  8231,  which 
is  a  175  W,  120  V,  60  Hz,  4-pole  machine.  The  system  parameters  were  used  to 
effectively  model  the  system  in  Simulink.  Both  the  stator  and  rotor  of  that  machine 
consist  of  3-phase  wye-connected  windings  and  a  stator  to  rotor  turns  ratio  Ns/Nr  =10. 

To  have  an  input  to  determine  the  values  of  the  estimated  parameters  required  use 
of  the  values  that  were  measured  in  [28]  after  a  DC  test,  a  no  load  test,  and  a  blocked 
rotor  test  on  the  generator  mentioned  above.  The  results  are  given  in  Table  2. 
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Table  2.  Measured  DFIG  parameters.  From  [28]. 


Parameters  0(t) 

Values 

rs 

12  [Cl] 

rr 

15[Q] 

Lls 

0.024 1[H] 

L,r 

0.024 1[H] 

Lm 

0.3342[H] 

2.  Simulation  Results 

After  conducting  a  large  number  of  simulations,  varying  the  conditions  and  some 
parameters  of  the  system,  it  turned  out  that  only  the  first  two  methods  gave  satisfactory 
results.  The  simulation  using  RLS-EF  method  was  leading  the  estimator  to  an  unstable 
output  after  a  period.  The  next  section  describes  the  results  of  the  first  two  methods. 

a.  Simulation  Using  LMS 

The  results  for  the  parameter  convergence  using  LMS  are  shown  in 
Figures  14,  15  and  16,  where  the  initial  condition  for  the  estimator  was  set  to  zero.  The 
sampling  time  for  all  the  simulations  was  set  equal  to  0. 1  ms. 
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*-m  0"0  n  (Ohms)  r$  (Ohms) 


Figure  14.  Convergence  of  rs  and  r,'  using  LMS  from  simulation. 


Figure  15.  Convergence  of  Lm using  LMS  from  simulation. 
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Figure  16.  Convergence  of  L/s  and  L /,.  using  LMS  from  simulation. 

One  can  see  that  the  parameter  estimates  converge  to  the  expected  values, 
but  it  takes  a  while  for  them  to  settle  down.  This  is  a  result  that  depends  on  the  non  zero 
values  of  the  adaptation  matrix  given  in  Table  1.  If  one  sets  the  values  of  initial  condition 
vector  closer  to  the  measured  values,  then  the  convergence  is  much  faster. 

b.  Simulation  Using  RLS 

The  results  of  the  simulation  using  RLS  are  shown  in  Figures  17,  18  and 
19. 
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Figure  17.  Convergence  of  rs  and  r/ using  RLS  from  simulation. 
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From  Figures  17  to  19,  one  can  conclude  that  the  parameter  convergence 
using  RLS  is  almost  instantaneous,  but  for  the  inductances,  the  estimates  oscillate  around 
the  expected  values.  On  the  other  hand,  the  estimates  for  the  resistances  are  extremely 
close  to  the  experimentally  measured  parameter  values. 

E.  CHAPTER  SUMMARY 

In  this  chapter,  different  algorithms  for  system  identification  were  presented 
analytically  and  simulated.  Based  on  the  results  and  given  the  fact  that  the  nature  of  the 
system  is  relatively  slow,  one  can  say  that  the  LMS  is  more  suitable  in  this  application 
since  it  performs  well  and  provides  relatively  fast  convergence.  In  the  next  chapter,  a 
basic  introduction  to  adaptive  control  theory  and  the  design  and  simulation  of  two 
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different  adaptive  controller  approaches,  along  with  simulation  results  under  different 
operational  conditions  of  the  system,  will  be  presented.  A  comparison  with  simulation 
results  obtained  from  classical  control  approaches  will  also  be  presented. 
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IV.  ADAPTIVE  CONTROL 


A.  INTRODUCTION 

According  to  Webster's  dictionary,  the  definition  of  adapt  is  "to  change  (oneself) 
so  that  one's  behavior  will  confonn  to  new  or  changed  circumstances"  [22].  According  to 
another  definition,  an  adaptive  system  is  any  physical  system  that  has  been  designed  with 
an  adaptive  viewpoint  [23].  From  the  above  expressions,  one  can  say  that  an  adaptive 
controller  is  a  controller  that  has  the  ability  to  change  its  behavior  to  adapt  to  changes  in 
the  systems  dynamics.  That  means  that  to  be  adaptive,  the  controller  has  to  have 
adjustable  parameters. 

Research  on  adaptive  control  systems  was  introduced  during  the  1950s,  mostly  for 
the  design  of  high  performance  flight  control  systems.  For  many  years,  there  has  been 
skepticism  among  scientists  about  the  use  of  adaptive  control,  since,  in  spite  of  the  fact 
that  the  adaptation  feature  sounds  very  attractive,  its  performance  and  reliability  may  be 
questionable.  In  particular,  when  one  applies  it  to  a  dynamic  system,  the  designer  should 
be  aware  that  adaptation  adds  another  loop,  which  might  make  the  overall  system 
unstable.  Some  important  issues  are  the  following  [29]: 

•  Because  of  internal  modulation  effects,  fast  adaptation  produces  a  high 
frequency  input.  This  input  excites  unmodeled  dynamics  easily  and  causes 
systematic  instability.  Typical  examples  from  [30]  have  shown  that  under 
very  mild  conditions  of  unmodeled  dynamics,  an  adaptive  control  system 
can  easily  become  unstable. 

•  There  is  no  guarantee  of  the  perfonnance  of  the  closed  loop  system  during 
the  transient  response. 

In  the  last  few  years,  new  applications  of  adaptive  control  have  started  to  appear, 
and  some  of  these  issues  have  been  addressed  with  a  degree  of  success  [31]. 
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An  adaptive  control  system  can  be  thought  as  having  two  loops.  One  loop  is  a 
normal  feedback  loop,  with  the  process  and  the  controller.  The  second  is  the  parameter 
adjustment  loop. 

The  system  is  modeled  in  terms  of  a  vector  of  parameters  9  that  may  be  known  or 
approximately  known.  The  goal  of  the  controller  is  to  continuously  produce  a  vector  of 
estimates  for  the  above  parameters  9  and  control  the  output  vector. 

B.  ADAPTIVE  CONTROL  THEORY 

There  are  many  different  approaches  to  adaptive  control.  The  model  reference 
adaptive  control  (MRAC)  and,  more  specifically,  the  direct  MRAC  (DMRAC),  is  used  in 
this  thesis.  In  this  control  scheme,  the  system  is  modeled  in  terms  of  a  vector  of 
parameters  0C  that  may  be  known  or  approximately  known.  This  parameter  vector  is 
updated  by  an  adaptive  law  using  input  and  output  data.  A  general  view  of  the  direct 
MRAC  is  shown  in  Figure  20. 


Figure  20.  Direct  MRAC  overview.  After  [31]. 


As  shown  in  Figure  20,  the  output  of  the  system  z(t)  is  compared  to  the  output  of  a 
reference  model  zref(t).  The  reference  model  is  chosen  to  generate  the  desired  trajectory 
for  the  plant  output  to  follow  [22].  The  controller  parameters  are  updated  continuously, 
depending  on  the  output  error: 
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z(t)  =  z(t)-zref(t).  (0.0) 

The  goal  of  the  controller  is  to  continuously  estimate  the  parameters  giving  a 
vector  6[t)  to  drive  the  error  of  the  output  vector  z(t)  to  zero.  The  rest  of  this  chapter 
examines  the  application  of  this  to  the  dynamic  model  of  the  DFIG. 

C.  ADAPTIVE  CONTROL  OF  DFIG 

Designing  a  control  scheme  required  the  evaluation  of  two  different  control 
approaches  according  to  which  part  of  the  system  should  be  adaptive: 

•  First  approach:  adaptive  control  of  torque  and  reactive  power  using  as 
output  the  rotor  currents  iqr  and  id,- 

•  Second  approach:  adaptive  control  of  torque  and  PI  control  of  the  reactive 
power  using  as  output  the  rotor  currents  iqr  and  id,- 

These  approaches  are  described  in  the  following  sections  and  the  results  obtained 
from  simulation  are  discussed. 

1.  First  Approach:  Adaptive  Control  for  Torque  and  Reactive  Power 

In  this  approach,  the  controller  was  designed  using  adaptive  control  for  both  rotor 
currents  components,  iqr  and  id,-,  which  form  the  z(t)  vector.  The  overall  design  is 
described  in  the  following  section. 

a.  Dynamic  Model 

The  voltage  equations  of  a  DFIG  expressed  in  an  arbitrary  chosen  qdO 
reference  frame  were  described  in  Chapter  II. 

For  the  purposes  of  this  design,  one  should  assume  that  the  reference 
frame  is  detennined  by  the  line  voltage  of  the  stator.  In  this  frame,  the  stator  voltage  is  a 
constant  with  value: 

vs{t)  =  vqs(t)-jvJt)  =  jVs  (0.0) 
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In  addition,  one  can  assume  the  stator  resistance  to  be  negligible;  i.e., 


rs  =0.  (0.0) 

Then,  defining  the  rotational  frequency  of  the  frame  co  to  be  the  same  as 
the  electrical  rotational  frequency  of  the  stator  coe,  one  can  write  Equation  (2.28)  as: 

pAs=-jcoAdqs+jVs.  (0.0) 


V 

The  solution  of  Equation  (4.4),  if  the  initial  condition  is  As(0)  =  — ,  is  a 

co 


constant  given  by 


V 

W)  =  ~- 

co 


(0.0) 


Since  the  flux  is  complex  As (t)  -  A .  (t)-y'/lds(t),  this  yields 


V 

A  (t)  =  —, constant 

co 

4(f)  =  o 


(0.0) 


The  fact  that  the  flux  is  a  constant  in  this  reference  frame  leads  to  the 
definition  of  a  constant  magnetizing  current  ims.  The  goal  of  the  rest  of  this  section  is  to 
rewrite  the  dynamic  model  of  the  DFIG  in  terms  of  the  rotor  current  components  iqr  and 
id,-  and  the  magnetizing  current  ims  only. 

Therefore,  using  Equations  (2.31)  from  Chapter  II  relating  the  fluxes  with 
the  currents,  one  can  obtain: 

K  _  +  LM  lm 

y  I  I  f  I  i 

_r  J  L  LM  Llr  “r  LM 


(0.0) 


Substituting  Equations  (4.6)  in  (4.7),  one  gets: 
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(0.0) 


lls  Lm  )iqs  +  LMiqr 


00 


(Lls+U'ds+U’dr  “° 


Solving  Equation  (4.8b)  for  ,  one  gets: 


l,  +  lm 


(0.0) 


Now  defining  the  constant  magnetizing  current  ims  as 

I  :  _  ; 

ms  qfs  'ms  ,  ‘ms  .  , 

lm  (>K 

and  substituting  Equation  (4.10)  in  Equation  (4.8a)  and  solving  for  iqs,  one  obtains: 


(0.0) 


ms 


V. 


co 


'  (Lis  +  LM)iqs  +  LMiqr 


(0.0) 


Lls  +  LM 


<ms~ 


Lls  +  LM 


qr 


Now  one  can  easily  rewrite  all  the  equations  in  tenns  of  ims  and  ir. 
Substituting  Equations  (4.9),  (4.10)  and  (4.11)  in  Equations  (4.7)  and  expanding  the 
expressions  for /l?rand/l one  obtains: 

\r  =  LMiqs  +  ((-/r  K qr  ^ 

lr=U: 


1m _ / - Id _ !'  )  +  (l  +l  )j' 

“ qr  -/M'.  .  'ms  .  .  1  qr  I  ~  'Hr  ~  UM  I1  qr 

‘-Is  +  ‘-M  ‘-Is  +  ‘-M 


>1-=-^— L,+(  1-  L 


L/s 


(Lls  +  LM)(Llr  +  Lm) 


)(4  +LM)'qr 


(0.0) 


and 


\ dr  LMids  +  (f/r  +  LM)idr 


2  =  I  (—  M 

dr  LM\  ,  ,  /'dr  1  V^/r  1  '-M/'dr 

Lls+LM 


-)idr+(Llr+LM)'d 


(0.0) 


^=(1-  f  T  “  T  Mr+LJdr 

(L+LM)(Llr+LJ 
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The  flux  linkages  equations  can  then  be  summarized  from  Equations  (4.6), 
(4.12)  and  (4.13)  as: 


4^ 

CO 

4=0 


4  -  .  M  .  'ms  +  °4r  +  LM  Vqr 

^-/S  +  *~M 

\ dr  ~  ^  +  LM)idr 


(0.0) 


where  a  is  defined  as  a  real  parameter  given  by  the  following  expression: 


<j  =  1-- 


(Lls  +LIVl)(Llr  +Lm ) 


,0  <  a  <  1 . 


(0.0) 


In  most  cases,  Lm  is  significantly  larger  than  L/v  and  Li,-,  which  means  that 
cr  is  very  small. 

From  Equations  (2.28)  in  Chapter  II,  after  ignoring  the  0  components, 
which  are  negligible  and  defining 


co  -  co  r=  co , 


(0.0) 


one  gets: 


P4=“fV“®4+V 

p4  =  -'4  +  ®4+'/*' 


(0.0) 


Substituting  Equations  (4.14)  in  Equations  (4.17),  one  finally  obtains  the 
desired  dynamic  model  of  DFIG  in  terms  of  magnetizing  and  the  rotor  currents: 


p/  = - ! - 


r.  1 

'qr-™dr+- 


°(Lr+Lb 


qr 


(0.0) 


P'dr  =  - 


V(Lr+L, 


-ldr+OJIqr+CO 


lms+' 


&(Lr  +  Lm  )(Ls  +  Lm  )  a(Lr  +  /.j 


'dr 


Now,  define  ao,ai,b  the  following  parameters: 
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(0.0) 


°o=- 


o,  = 


^r+U 

l  2 


a(Lr  +Lm)(Ls+LI IV 
1 


-'ms  ' 


b  = 


°(Lr+L, 


Substituting  Equations  (4.19)  in  Equations  (4.18),  one  ends  up  with  two 
simple  first  order  models: 

Piqr=(]0iqr-o)i'dr+bvqr 

P'.=°o'>®V+®Oi+h' 


(0.0) 


By  defining  the  state  vector  z  = 


V 

and  the  control  vector  v  = 

V 

J  dr  _ 

Vdr_ 

and 


choosing  two  arbitrary  positive  constants  am,  [>m,  one  can  rewrite  Equations  (4.20)  in 
terms  of  the  state  vector  and  a  desired  dynamic  model: 


l(t)  =  -am  z(t)  +  j3mv(t)  +  ®l(t)0  +  a) 


-I 


dr 


qr 


(0.0) 


with 


follows 


0T(t)  = 


0  = 


V  0  V 

• dr  ®  Vdr 


«0  +am 

a1 

b-Pm 


(0.0) 


The  significance  of  Equation  (4.22)  is  that  by  defining  the  control  signal  as 


'dr 


-OT(t)0(t)  +  Wref(t) 


(0.0) 
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with  9[t )  an  estimate  of  the  system  parameter  vector  0(t) ,  the  closed  loop  system 
becomes 

z(t)  =  - amz(t )  +  <bT(t)£(t)  +  Wref(t)  (0.0) 

where 

d(t)  =  0-0(t)  (0.0) 

In  the  ideal  case  of  perfect  knowledge  of  the  DFIG  dynamics,  the 
parameter  error  vector  0{t)  in  Equation  (4.25)  is  zero,  and  the  state  vector  z(t)  tracks  the 
reference  model: 

Iref  (f)  =  -am  Zref  (t)  +  Wref  (t)  (0.0) 

The  tracking  error  z(t)  is  defined  as: 

z(t)  =  z(t)-zref(t).  (0.0) 

Then,  combining  Equations  (4.24)  and  (4.26)  and  substituting  in  Equation 
(4.27),  one  ends  up  with  the  following  expression: 

z(t)  =  -am~z(t)  +  0Tmt).  (0.0) 

Now,  estimating  the  controller  parameters  requires  the  development  of  a 
differential  equation,  called  the  update  law.  That  differential  equation  is  of  the  form: 

9(t)  =  F{y(t),  e(t),  9(t),  O(t)} .  (0.0) 

Using  the  steepest  descent  algorithm,  we  ensure  that  the  parameter  error 

decreases  monotonically  with  time.  In  particular  the  definition  of  0(t)  is  as  follows: 

6[t)  =  //fi(t)z(t) ,  (0.0) 

where  ju>0  is  the  adaptation  gain  of  the  controller.  The  definition  of  the  following 
Lyapunov  function  is: 


46 


v[e,t)=\  ~z[t)2 +^~  d(t)2.  (o.o) 

2  2  JU 

As  also  stated  in  Chapter  III,  for  the  system  to  be  stable,  the  Lyapunov 
function  described  in  Equation  (4.31)  has  to  be  positive  definite,  and  its  derivative  with 
respect  to  time  must  be  negative  semi-definite.  By  definition,  it  is  positive  definite,  since 
in  this  case  the  parameters  are  positive  real  numbers  and  the  norm  of  the  parameter 
vectors  are  used.  The  derivative  with  respect  to  time  of  V(0(t))  is 

vm))=zT(t)z(t)+-0(tm),  (o.o) 

a 

where,  since  6{t)  ~  0 ,  one  can  use  the  fact  that: 

m=-m.  (o.o) 

In  Equation  (4.32),  if  one  substitutes  Equations  (4.28)  and  (4.30),  one  gets: 

mt) = zJ(t)(-aml(t) + ®Tmt))-mt)~z(t)  ym = -«m  ||i(t)||2 .  (o.o) 

From  Equation  (4.34),  V(6, t)  is  negative  semi-definite  for  every  value  of 

time  t. 

Equation  (4.34)  implies  that: 

V(t)  <  V(0 )Vt 


1  ~  2  1  ~  2  1  ~  2  1  ~  2 

-m  +—m  <-Z( o)  +—^(o) 

2  2ju  2  2ju 


which  means,  that  if  one  chooses 

z(0)  =  z(0)^z(0)  =  0, 


And  from  Equation  (4.35) 


(0.0) 
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(0.0) 


z(t) 


< 


m 

n/a 


Equation  (4.37)  means  that  z(t)  is  bounded. 


b.  Simulation  Model 

A  Simulink  model  for  the  controller  was  created  using  the  above  equations 
for  the  system.  The  inputs  to  the  controller  are  the  following  quantities,  all  expressed  in 
the  abc  reference  frame:  stator  voltages  Va/,cs,  stator  and  rotor  currents  iabcs  and  iabcr 
respectively,  cor  and  the  reference  currents  iqrref  and  i dr  ref  The  outputs  of  the  controller 
are  the  rotor  voltages  that  go  into  the  DFIG  model.  A  general  block  diagram  of  the 
controller  is  shown  in  Figure  2 1 . 


Reference  Model 


Theta_hat(t) 

Vref(t) 

Vqdr_control 
omega_tildl-idr  iqr] 

Phi_trans(t) 

CONTROL 


Vqdr_control 


Figure  21.  General  view  of  the  first  controller  in  Simulink. 
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This  model  was  created  using  Equations  (4.18)  to  (4.37).  A  detailed  view 

of  the  controller  block,  which  gives  the  values  of  the  control  vector  [  according 

to  Equation  (4.23),  is  shown  in  Figure  22.  The  output  is  then  transformed  back  to  the  abc 
reference  frame  to  supply  the  rotor  with  the  proper  control  voltage. 


Figure  22.  Control  block  for  the  first  approach. 

The  reference  model  shown  in  Figure  21  has  a  first  order  transfer  function 

of  the  form 

T(s)  =  — .  (0.0) 

s  +  °m 

The  control  goal  is  to  independently  control  the  rotor  currents  to  minimize 
the  reactive  power  of  the  DFIG  and  to  produce  the  maximum  real  power  for  every 
different  value  of  the  wind  speed  respectively. 

From  Equation  (2.36),  after  substituting  the  values  of  Vqs  and  V(/s  from 
Equation  (4.2),  one  can  obtain  the  following  expression  for  the  reactive  power: 

(o.o) 

which  means  that  Qe  is  proportional  to  iqs,  so  to  have  Qe  =  0  ,  iqs  has  to  be  equal  to  zero; 
i.e., 
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(0.0) 


qs_ref 


=  0  . 


From  equation  (4.11),  substituting  the  above  value  for  iqs  ref  from 
Equation  (4.40),  one  gets: 


qr_ref  ^ ms  ? 


(0.0) 


which  means  that  iqr  should  be  equal  to  the  magnetizing  current,  which  is  constant  after 
the  system  reaches  steady-state,  and  is  given  by  Equation  (4.10). 


On  the  other  hand,  keeping  the  power  output  maximum  requires  that  the 
electric  torque  of  the  DFIG  be  equal  to  the  aerodynamic  torque  that  gives  the  maximum 
power  output  of  the  WT  for  each  wind  speed,  which  is  given  by  [12] 


where 


'  e-ref 


_  t  —  is  o 

aero  _optimal  'opt^^r 


(0.0) 


opt 


(0.0) 


and  Q,  is  the  rotational  speed  of  the  WT  in  revolutions  per  minute  (RPM).  This  is  given 
by 


Q  = 


Gear  Ratio  P 


-  co. 


(0.0) 


where  Gear  Ratio  is  the  gear  box  ratio,  P  is  the  number  of  poles  of  the  DFIG,  and  cor  is 
the  electrical  rotational  speed  of  the  rotor  in  radians  per  second. 

The  expression  of  the  electric  torque  output  of  the  DFIG  is  given  by 
Equation  (2.35).  Substituting  the  values  of  the  flux  linkages  from  Equation  (4.6),  one  gets 


T=- 


( 31 

f 

UJ 

UJ 

u 


v. 


CO „ 


•dsi’ 


(0.0) 


and  by  substituting  z*  from  Equation  (4.9)  into  Equation  (4.45),  one  obtains 
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(0.0) 


T=- 


'3' 

' P ) 

1/ 

K 

|  TVJ 

v2v 

v2j 

\lm 

6)b 

^~ts  l~M 

dr  ' 


The  goal  is  to  track  a  reference  torque,  which  requires  Te  =  Te  ref .  That 

means  that  from  Equations  (4.44)  and  (4.46)  one  should  set  the  reference  current  id,-  that 
is  driving  the  reference  model  to  be 


'  e_ref 


{  3l 

UJ 

UJ 

V ■ 


(0.0) 


Running  the  model  required  the  inclusion  of  a  block  that  takes  the  current 
wind  speed  as  input  and  gives  the  instant  aerodynamic  torque  as  output,  which  is  the 
excitation  of  the  DFIG.  The  WT  turbine  speed  is  then  given  by 

A=Taero-Te,  (0.0) 

where  J  is  the  combined  inertia  of  the  WT  and  the  DFIG  rotor.  In  this  thesis,  the  rotor 
inertia  was  neglected  due  to  its  very  small  value  compared  to  the  inertia  of  the  WT. 

The  comparison  of  the  control  approaches  presented  here  with  other 
controllers  required  the  use  of  the  parameters  of  the  DFIG  and  the  WT  described  in  [7], 
which  are  shown  in  Table  3. 


Table  3.  WT  and  DFIG  parameters.  From  [7], 


Wind  Turbine 

DFIG 

R  =3.8  m 

Pmax= 50  hp= 37.285  kW 

J= 3.362  kgm2 

P= 4 

a=19.346 

Rs= 0.082 

6=9.41 17 

Rr= 0.228 

o 

CM 

II 

O 

Ts=0.0355  Hy 

Gear  Ratio  16:1 

Tr=0.0355  Hy 

/loPt=6.4 

Zm=0.0347  Hy 

cob=_= 2x60  rad/s 

Vs  =380  V 
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c.  Simulation  Results 

After  a  large  number  of  simulations,  the  system  parameters  that  gave  the 
best  results  are  shown  in  Table  4.  The  sampling  time  for  the  simulation  was  chosen  equal 
to  0.8  ms.  For  these  values,  the  system  was  stable  and  tracked  the  reference  model  for 
wind  speeds  from  3.2  m/s  to  10.8  m/s,  which  is  a  good  range.  This  covered  almost  all  the 
operational  range  of  the  WT. 


Table  4.  Final  simulation  parameter  values  for  first  approach. 


Parameter 

Value 

100 

bm 

250 

fill 

0.75 

1-122 

50 

fi33 

0.4 

In  Table  4,  it  can  be  seen  that  the  adaptation  gain  for  the  parameter 
estimation  was  chosen  to  be  a  diagonal  matrix  as  seen  in  Chapter  III  to  improve  the 
convergence  of  the  estimates.  The  initial  condition  for  the  estimates  was  chosen  equal  to 
70%  of  the  expected  values  calculated  using  Equation  (4.22). 

The  wind  speed  profile  used  is  shown  in  Figure  23  along  with  a  plot  of  the 
absolute  values  of  the  actual  output  torque  Te  and  the  reference  torque  Tref-  The  wind 
speed  profile  was  chosen  in  order  to  prove  the  performance  of  the  controller  even  for 
large  rapid  changes  in  wind  speed,  greater  than  5  m/s,  which  can  be  considered  wind 
gusts  according  to  [33]. 

The  output  torque  is  tracking  the  reference  value  for  a  range  of  wind 
speeds  from  3.2  m/s  to  10.7  m/s  as  shown  in  Figure  23. 
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Figure  23.  Wind  speed  and  torque  during  simulation  of  the  first  approach. 


Zooming  in  on  the  torque  plot  from  time  90  s  to  120  s  in  Figure  24,  one 
can  see  that  even  for  wind  changes  of  6.8  m/s,  the  difference  between  reference  and 
actual  output  is  less  than  1.65  %.  For  lower  wind  speeds,  the  difference  is  always  less 
than  0.7  %,  which  is  a  very  satisfactory  result. 


Figure  24.  Magnified  torque  plot  during  simulation  of  the  first  approach. 


The  plots  of  the  real  power  Pe  and  reactive  power  Qe  during  the  simulation  are 
shown  in  Figures  25  and  26,  respectively. 
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Figure  25.  Real  power  (Pe)  during  simulation  of  the  first  approach. 


Figure  26.  Reactive  power  ( Qe )  during  simulation  of  the  first  approach. 


In  Figures  25  and  26,  one  can  see  that  the  ratio  of  Qe  over  Pe  is  very  small 
and  remains  less  than  0.2%  for  most  the  simulation  time.  This  result  indicates  that  the 
second  goal  of  the  design  is  also  satisfied.  The  only  time  that  Qe  increases,  but  still  with 
values  less  than  1%  of  Pe,  is  during  extreme  changes  in  wind  speed,  and  that  happens  for 
a  very  short  period. 

The  rotor  currents  during  the  simulation  are  shown  in  Figure  27. 
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Figure  27.  Comparison  of  the  rotor  currents  iqr  and  /)/,■  with  their  reference  values 

during  simulation  of  the  first  approach. 


2.  Second  Approach:  Adaptive  Control  of  Torque,  PI  Control  of  the 
Reactive  Power 


A  second  approach  was  attempted  to  simplify  the  controller  with  only  one 
adaptive  loop  for  torque  control.  The  difference  in  this  approach  is  that  instead  of 
controlling  both  iqs  and  id,-  adaptively,  a  PI  controller  was  used  for  iqs  and  the  adaptive 
control  for  id,-  was  kept. 


a.  Dynamic  Model  of  the  Controller 


Equations  (4.2)  to  (4.18)  are  valid  for  this  approach.  From  Equation 
(4.1 1),  solving  for  iqr,  one  gets 


/  =i  -LM+Llsj 
qr  ms  .  qs 

‘-M 


(0.0) 


Substituting  Equation  (4.49)  for  the  last  one  of  Equations  (4.18), 

one  obtains 
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1 


P'dr=- 


rLrr'*+^-^J  +  ® 


+/-m) 


'ms  +■ 


CT(/-/r  +  LM  )  (*-b  +  LM  )  Cr(ifr  +  LM  ) 


.  (0.0) 


Pi*=~  7~  T  -/dr-^(^L)/  +<5(i  — 
oK+U  LM 


a(Llr+LM)(Lls+LJ  ms  <t(4+Lm)  * 


One  defines  the  following  parameters: 


°o=- 


°K+iJ 


°1  = 


Lk  +  Lm 


a  =  (1 - ^ - )/ 

a^+ZJ^+U 

i 

b  =  - 


°(Llr+U 

Substituting  Equations  (4.51)  in  Equation  (4.50),  one  gets 
Pi dr  =  °o idr  ~  °i  )  +  o2®  ■ +  bvdr . 


(0.0) 


(0.0) 


Now  zr  =  idr  and  y  =  i/dr  are  called  the  state  and  control  signal, 

respectively,  which  in  this  approach  are  scalars.  Then,  for  any  arbitrary  constant  am,  fi„ 
one  can  write  Equation  (4.52)  in  the  form 

z(t)  =  -amz(t)  +  /3mv(t)  +  ®T(t)0.  (0.0) 

where 

®T(t)  =  [idr  ~coiqs  a  vdr]  (0.0) 

and 


0 


«0 

°2 

b~Pm 


(0.0) 


In  this  case,  the  control  signal  is 
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/?mv(f)  =  -®T(t)£+!/!/re/(f)- 


(0.0) 


Now  substitute  Equation  (4.56)  into  Equation  (4.53),  and  one  gets: 

z(t)  =  -amz(t)  +  ® T(t)6(t)  +  Wref(t) ,  (0.0) 

where  0[t)  is  given  by  Equation  (4.25). 

b.  Simulation  Model 

The  general  view  of  the  controller  in  Simulink  is  shown  in  Figure  28. 
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Figure  28.  General  view  of  the  second  controller  model. 


Equations  (4.30)  to  (4.37),  which  were  used  in  the  first  approach  for  the 
parameter  estimates,  were  also  used  for  this  case. 

For  iqs,  a  PI  controller  was  used  with  the  transfer  function 

K^s  +  K^.  (0.0) 

s 
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where  Kq,  K\  are  the  gains  of  the  PI  controller.  The  input  of  this  controller,  the  error 
between  the  reference  and  the  actual  measured  value  of  iqs  and  the  output  is  vqr. 

Furthermore,  Equations  (4.42)  to  (4.48)  that  give  the  expression  for  Ur  ref 
and  Q,  were  again  used  here. 

In  this  approach,  the  output  of  the  adaptive  controller  is  vdr  according  to 
Equation  (4.57);  two  low  pass  (LP)  fdters  were  added  in  the  output  of  the  adaptive  part  to 
filter  higher  frequency  components  that  might  affect  the  stability  of  the  system.  A 
detailed  view  of  the  control  block  diagram  is  shown  in  Figure  29. 


Figure  29.  Control  block  for  the  second  approach. 

The  cutoff  frequency  ac  in  the  LP  filter  shown  in  Figure  29  was  100,  and 
for  the  second  filter,  which  is  a  12th  order  Butterworth  filter,  the  cutoff  frequency  was  10 
times  higher. 


c.  Simulation  Results 

After  a  large  number  of  simulations,  the  system  parameters  that  gave  the 
best  results  for  this  approach  are  those  shown  in  Table  5.  The  sampling  time  for  the 
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simulation  was  equal  to  0.8  ms.  For  these  values,  the  system  was  stable  and  tracked  the 
reference  model  for  wind  speeds  from  3.2  m/s  to  10.8  m/s. 


Table  5.  Final  simulation  parameter  values  for  the  second  approach. 


Parameter 

Value 

a  in 

100 

bm 

250 

fill 

100 

In  Table  5,  it  can  be  seen  that  the  adaptation  gain  for  the  parameter 
estimation  in  this  case  was  a  single  value,  which  is  about  two  orders  of  magnitude  larger 
than  the  adaptation  gains  for  the  first  approach  given  in  Table  4.  This  value  was  chosen 
because  comparable  results  could  not  be  achieved  between  the  two  control  approaches 
with  smaller  values  of  //.  The  initial  condition  for  the  estimates  was  chosen  equal  to  70% 
of  the  expected  values  calculated  using  Equations  (4.56). 

The  wind  speed  profile  is  the  same  as  the  one  used  in  the  previous 
approach,  which  allows  for  a  comparison  between  approaches  in  the  next  chapter.  This 
plot,  along  with  a  plot  of  the  absolute  actual  output  torque  Te  and  the  reference  (optimal) 
torque  Tref,  is  shown  in  Figure  30. 
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From  Figure  30,  it  can  be  seen  that  the  output  torque  tracks  the  reference 
value  pretty  well  for  a  range  of  wind  speeds  from  3.2  m/s  to  10.7  m/s,  but  during  the 
sudden  wind  changes,  the  actual  output  torque  spikes  for  a  few  moments.  A  zoomed  in 
view  of  the  torque  plot  from  time  97  s  to  1 19  s  is  illustrated  in  Figure  3 1 . 


Figure  3 1 .  Magnified  torque  during  simulation  of  the  second  approach. 


It  can  be  seen  in  Figure  31  that  for  large  and  fast  wind  changes  that 
simulate  wind  gusts,  the  torque  tends  to  lose  tracking  for  few  moments.  This  is  especially 
true  when  the  wind  decreases,  e.g.  at  100  s  or  at  117  s  of  simulation  time,  where  the 
difference  between  the  two  values  even  reaches  100%  but  recovers  quickly. 

The  plots  of  the  real  power  Pe  and  reactive  power  Qe  during  the  simulation 
of  the  second  approach  are  shown  in  Figures  32  and  33,  respectively. 


Figure  32.  Real  power  (Pe)  during  simulation  of  the  second  approach. 
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Figure  33.  Reactive  power  ( Qe )  during  simulation  of  the  second  approach. 


It  is  can  be  seen  in  Figure  33  that  Qe  is  quite  high  compared  to  Pe.  This  is 
clearer  if  one  zooms  in  on  a  part  of  the  simulation,  as  shown  in  Figure  34. 


Figure  34.  Comparison  between  real  and  reactive  power  Pe  and  Qe,  respectively,  for 

part  of  the  simulation  during  the  second  approach. 


Finally,  the  plots  of  the  currents  iqs  and  ijr  are  shown  in  Figure  35. 
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Figure  35.  Comparison  of  the  currents  iqs  and  id,  with  their  reference  values  during 

simulation  of  the  second  approach. 

It  can  be  seen  in  Figure  35  that  id,-  is  tracking  its  reference  value  very  well, 
but  iqs  seems  to  lose  tracking  during  sudden  wind  speed  changes. 

3.  Performance  Comparison 

To  be  able  to  reach  a  conclusion  on  the  performance  of  the  controllers  designed  in 
this  study,  a  Simulink  model  of  a  classical  approach  for  the  control  of  a  WECS  was 
created.  This  is  described  in  [12],  again  using  the  parameters  of  the  system  given  in  [7]. 
The  model  was  created  for  [33]. 

Two  different  fault  cases  were  simulated  to  compare  the  robustness  of  the  control 
algorithms.  These  two  are: 

•  a  wind  gust,  i.e.,  a  sudden  increase  in  the  wind  speed  from  6  to  10  m/sec, 

•  doubling  of  the  value  of  the  rotor  resistance  r,.  to  demonstrate  what 
happens  during  a  parameter  variation. 
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The  wind  profile  shown  in  Figure  36  was  assumed. 


Figure  36.  Wind  profile  for  comparison  of  the  control  approaches. 


The  output  torque  of  the  three  different  control  approaches  during  simulation  is 
shown  in  Figure  37. 


Figure  37.  Comparison  of  output  torque  Te  and  reference  torque  Tref  during  the 

simulation  of  a  wind  gust. 


Magnifying  Figure  37  for  about  10  s  during  the  top  of  the  wind  gust  provides 
Figure  38. 
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Figure  38.  Magnified  plot  of  output  torque  Te  and  reference  torque  Tref  during  the 

simulation  of  a  wind  gust. 


From  Figure  38,  one  can  see  that  even  though  the  performance  of  all  the 
approaches  considered  is  very  good,  the  output  torque  of  the  first  approach  and  the 
classical  approach  with  the  PI  controllers  are  identical  and  almost  identical,  respectively, 
to  the  reference  torque  throughout  the  simulation,  while  the  second  approach  shows  some 
oscillation  during  sudden  wind  speed  changes.  Furthermore,  the  torque  of  the  classical 
approach  is  slightly  closer  to  the  reference  value,  but  has  more  overshoot  compared  to  the 
first  approach. 

The  resistance  of  the  rotor  is  a  parameter  that  affects  the  shape  of  the  torque 
output  of  the  DFIG  [19]  and  is  illustrated  in  Figure  39. 
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Figure  39.  Typical  DFIG  torque  versus  time  curve  for  different  values  of  rr. 


For  the  purposes  of  this  thesis,  in  order  to  investigate  the  performance  of  the 
controllers  that  we  proposed  in  previous  sections,  the  rotor’s  resistance  rr  was  increased 
to  twice  its  initial  value  at  t=  1 5  s  of  the  simulation,  and  all  the  control  approaches  where 
simulated.  The  output  torque  plots  during  the  above  change  in  rr  from  the  simulations  are 
shown  in  Figure  40. 


Figure  40.  Comparison  of  the  variation  of  the  output  torque  curves  during  simulation 

of  a  variation  of  rr. 
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As  seen  in  Figure  40,  Te  for  the  second  approach  has  a  high  overshoot  and  even 
after  three  seconds  does  not  completely  track  the  reference  torque  Trej.  Between  the  other 
two  cases,  the  first  approach  appears  to  have  the  best  performance  since  it  has  a  small 
spike  in  Te  but  recovers  faster  than  the  PID  approach.  This  can  also  be  seen  in  Figure  41, 
where  the  speed  of  the  DFIG  for  the  three  approaches  is  presented  around  the  step  change 
in  rr;  it  is  clear  that  the  first  adaptive  approach  gives  the  smoother  speed  curve  among  the 
three  different  cases. 


Figure  41.  Comparison  of  the  variation  of  the  DFIG  speed  during  the  simulation  of  a 

variation  of  rK 


D.  CHAPTER  SUMMARY 

In  this  chapter,  the  analytical  design,  modeling  and  simulation  of  two  MRAC 
approaches  for  a  DFIG-based  WECS  were  presented.  The  simulation  results  were 
compared  with  results  from  simulation  of  a  more  classical  control  approach  for  this  kind 
of  system,  which  was  presented  in  [12].  The  conclusions  from  these  comparisons  along 
with  recommendations  for  future  research  are  discussed  in  Chapter  V. 
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V.  CONCLUSIONS  AND  RECOMMENDATIONS 


A.  SUMMARY  -  CONCLUSIONS 

The  problem  of  the  control  of  a  doubly-fed  induction  generator  based  wind  energy 
conversion  system  with  fixed  pitch  horizontal  axis  wind  turbine  was  addressed  in  this 
thesis.  The  adaptive  control  used  was  direct  model  reference  adaptive  control.  This 
research  effort  addressed  issues  relative  to  the  identification  of  the  system’s  parameters 
and  the  design  of  adaptive  control  techniques  in  order  to  compensate  for  uncertain  or 
time -varying  dynamics. 

In  particular  three  system  identification  algorithms  were  presented  based  on 
standard  techniques  such  as  least  mean  squares  (LMS)  or  recursive  least  squares  (RLS) 
with  and  without  exponential  forgetting.  The  significance  of  these  techniques  is  that  a 
dynamic  system  of  the  DFIG  can  be  properly  parameterized  on  the  basis  of  experimental 
results. 

Furthermore,  two  different  adaptive  control  approaches  were  introduced  and 
simulated  for  application  in  wind  power  generation.  In  this  case  it  turns  out  that  the 
parameterized  dynamic  model,  based  on  the  magnetic  fluxes  of  stator  and  rotor,  can  be 
simplified  considerably  and  only  a  few  parameters  have  to  be  estimated.  This  is  due  to 
the  stator  being  connected  to  the  grid  so  that  the  stator  voltage  is  constant  in  magnitude, 
frequency  and  phase.  Given  this  characteristic,  the  magnetizing  current  becomes  one  of 
the  parameters  to  be  estimated,  which  results  in  considerable  simplification  of  the 
dynamic  model. 

These  approaches  have  been  tested  in  numerous  simulations  and  seem  to  present 
an  attractive  solution  to  the  problem  of  compensating  for  uncertain  dynamics  such  as 
varying  resistances  and  impedances  in  the  circuitry. 
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B. 


REC OMMEND ATION S-FUTURE  RESEARCH 


Extensive  computer  simulations  demonstrate  that  the  algorithms  presented  for 
adaptive  control  work  as  expected.  Next  step  in  this  research  will  be  the  implementation 
of  the  proposed  controllers  in  a  laboratory  environment  using  a  direct  current  motor  to 
simulate  the  wind  turbine  and  building  the  controllers  using  Xilinx  software  in  Simulink 
and  field  programmable  gate  arrays  to  establish  bidirectional  communication  between  the 
software  and  the  equipment. 

This  effort  will  involve  not  only  actual  real-time  implementation  but  also  the 
design  of  an  appropriate  model  for  the  wind  profile  so  that  wind  effects  can  be  simulated 
by  the  generation  of  appropriate  torque  by  the  DC  motor.  The  ultimate  benefit  of  this 
research  will  be  more  efficient  wind  power  generators. 
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APPENDIX:  MATLAB  CODES  -  SIMULINK  MODEL 


In  this  Appendix,  the  code  that  was  used  for  the  initialization  of  the  simulations, 
as  well  as  the  basic  block  of  the  Simulink  models  for  the  controllers  are  presented. 

A.  MATLAB  CODE  FOR  INITIALIZATION  OF  THE  SIMULATIONS 

%%  INITIAL  CONDITIONS  FILE 
clc; 

clear  all; 

omega_b  =  2*pi*60; 
tstep= . 0005 ; 
tstop=120; 
twopiby3  =  2*pi/3; 
poles  =  4; 

J=3.362;  %  kg  mA2 

polesby2J  =  poles/2/J; 

GearRatio=16; 

%Parameters  from  paper  on  nonlinear  control 

rrO  =0.228;  %  ohms 

rs=0.02;  %  ohms 

Xls=omega_b* ( 0 . 0355-0 . 0347 ) ;  %  ohms 
Xm  =omega  b* 0.034 7;  %  ohms 

Xlr  =omega  b* (0 . 0355-0 . 0347) ;  %  ohms 

Vs=380;  %  Volts  Grid  phase  voltage 
omega  e=omega  b;  %electrical  frequency 

%Parameters  -  Equivalent  circuit 

Ls=Xls/ omega_b; 

Lm=Xm/ omega  b; 

Lr=Xlr/omega  b; 
rsbyXls  =  rs/Xls; 
rrbyXlr  =  rrO/Xlr; 

Xaq  =  1/ ( 1/Xm+l/Xls+l/Xlr) ; 

Xad  =  Xaq; 

XaqbyXls  =  Xaq/Xls; 

XaqbyXlr  =  Xaq/Xlr; 

XadbyXls  =  Xad/Xls; 

XadbyXlr  =  Xad/Xlr; 

V_phase  =  380*sqrt (2) /sqrt (3) ;  %Peak  value  used  in  Simulink  for  source 
model 
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start  at  synchronous  speed 


%Steady  state  analysis 
omegar  ic  =  omega  b*.85; 
slip  =  0.001; 


Z1  =  rs  +  j*Xls; 

Z2  =  j  *Xm; 

Z3  =  rrO/slip  +  j*Xlr; 

las  =  V  phase/sqrt (2) / (Z1  +  Z2*Z3/ (Z2  +  Z3) ) ; 

Iar  =  -las* (Z2*Z3/ (Z2+Z3) ) /Z3; 

Te  =  3*poles/2*Xm/omega_b*real ( j *con j ( las ) * Iar ) ; 
Iqse  =  sqrt (2) *real (las) ; 

Idse  =  -sqrt (2) *imag (las) ; 

Iqre  =  sqrt (2 ) *real ( Iar) ; 

Idre  =  -sqrt (2 ) *imag ( Iar)  ; 
psi  mqe  =  Xm* ( Iqse+Iqre) ; 
psi  mde  =  Xm* ( Idse+Idre) ; 
psi  qse  =  Iqse*Xls  +  psi  mqe; 

psi  dse  =  Idse*Xls  +  psi  mde; 

psi  qre  =  Iqre*Xlr  +  psi  mqe; 

psi  dre  =  Idre*Xlr  +  psi  mde; 

psi_qsic=psi_qse; 
psi_dsic=psi_dse; 
psi_qric=psi_qre; 
psi  dric=psi  dre; 


%WIND 

rho=1.225;  %  air  density 

lam=0 . 1 : .1:10; 

C= (19. 346. /lam) .* (9. 4 117. /lam  -1) . *exp (-20 . /lam) 
%  plot (lam,  C) 
lam_opt=6 . 4 ; 

C_opt=C  (  64 )  ; 

K_opt=rho*pi* (3 . 8A5) *C_opt/ (2*lam_optA2) ; 
lamda=0 . 8 ; 


%MRAC  parameters 
am=100.0;  %  this  works 

bm=250; 

%  Filter  frequency 
ac=500 ; 

sigma=l- (LmA2) / ( (Ls+Lm) * (Lr+Lm) ) ; 
b=l/ (sigma* (Lr+Lm) ) ; 
a0=-rr/ (sigma* (Lr+Lm) ) ; 
c=50  ; 

%  controller  parameters 

im0=Vs/ (omega  e*Lm) ;  %  estimated  magnetizing 

%  true  values 

theta ( 1 ) =a0+am; 
theta (2 ) =Lm/ (Ls+Lm) ; 
theta (3) =-im0*a0*Lm/ (Ls+Lm) ; 
theta ( 4 ) =b*Lm/ (Ls+Lm)  -  bm; 
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theta (5)  =  (Ls+Lm) /Lm; 

theta (6) = (1+ (LmA2) / (sigma* (Ls+Lm) * (Lr+Lm) ) ) *imO; 
theta ( 7 ) =b-bm; 

%  Estimation  bounds 

rho  theta=50;  %  parameter  error  in  "percent" 
d_theta=abs (theta) *rho_theta/100; 
theta  min=theta-d  theta; 
theta  max=theta+d  theta; 

%  initial  estimates 

thetaO=theta+d_theta . * (2*rand (1, 7) -1) ; 

B.  SIMULINK  BLOCKS 


1.  First  Adaptive  Control  Approach 


a.  Overall  System  View 


Grid 


WINO  SPEED  PROFILE 


■vT 


Aerodynamic 

Torque 


Oes 


vcs 


r*  Ta 


DFIG 


Pes 


To  Worts  Dace  1 


lambda 
(opt  =  6.4) 


Gsi- 

Conversion  of  rad/s 
toRPM 
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b.  Aerodynamic  Torque  Input 


UA2 


c.  General  View  of  the  Controller 
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d.  Preprocess 


e.  Parameter  Estimation 
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f  Control 


2.  Second  Approach 

The  above  Simulink  blocks,  in  sections  a,  b,  c,  e,  are  common  for  both 
approaches.  Those  which  are  different  for  the  second  approach  follow. 

a.  Preprocess 
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b. 


Control 


75 


THIS  PAGE  INTENTIONALLY  LEFT  BLANK 


76 


LIST  OF  REFERENCES 


[1]  U.S.  Energy  Information  Administration  (EIA), ’’International  Energy  Outlook 
2010”  [Online].  Available:  http ://www.eia . gov/oiaf/ieo/ electricity.html  (accessed 
August  5,  2010). 

[2]  F.  Birol,  “Power  to  the  People,  The  World  Outlook  for  Electricity  Investment,” 
International  Atomic  Energy  Agency  Bulletin,  Bulletin  46/1,  pp.  9-12,  June  2004. 

[3]  REN21,  “Renewables  2011  Global  Status  Report,”  Technical  Report,  REN21 
Secretariat,  Paris,  2011. 

[4]  U.S.  Department  of  Energy,  Office  of  Energy  Efficiency  and  Renewable  Energy, 
“20%  Wind  Energy  by  2030,”  Technical  Report  No.  DOE/GO-102008-2567,  U.S. 
DOE,  Washington,  DC,  2008. 

[5]  K.  E.  Johnson,  “Adaptive  Torque  Control  of  Variable  Speed  Wind  Turbines,” 
Technical  Report  No  NREL/TP -500-36265 ,  National  Renewable  Energy 
Laboratory,  Golden,  Colorado,  August  2004. 

[6]  J.  I.  Jang,  Y.  S.  Kim,  and  D.  C.  Lee,  “Active  and  Reactive  Power  Control  of 
DFIG  for  Wind  Energy  Conversion  Under  Unbalanced  Grid  Voltage,”  Presented 
at  Power  Electronics  and  Motion  Control  Conference,  2006.  IPEMC  2006. 
CES/IEEE  5th  International,  Shanghai,  August  2006. 

[7]  F.  Valenciaga,  P.  F.  Puleston,  and  S.  K.  Spurgeon,  “A  Geometric  Approach  for 
the  Design  of  MIMO  Sliding  Controllers,  Application  to  a  Wind-driven  Doubly 
Fed  Induction  Generator,”  Int.  J.  of  Robust  and  Nonlinear  Control,  vol.  19,  issue 
1,  pp.  22-39,  January  2009. 

[8]  V.  Utkin,  J.  Guldner,  and  Shi  J.,  Sliding  Mode  Control  Theory  in 
Electromechanical  Systems,  Taylor  &  Francis,  London,  1999. 

[9]  A.  Sabanovic,  L.  Fridman,  and  S.  K.  Spurgeo,  “ Variable  Structure  Systems:  from 
Principles  to  Implementation IEEE  Book  Series.  IEEE,  London,  2004. 

[10]  S.  A.  Frost,  M.  J.  Balas,  and  A.  D.  Wright,  “Direct  Adaptive  Control  of  a  Utility- 
scale  Wind  Turbine  for  Speed  Regulation,”  Int.  J.  of  Robust  and  Nonlinear 
Control,  vol.  19,  issue  1,  pp.  59-71,  January  2009. 

[11]  A.  Martinez,  F.  Lescher,  and  J.  Y.  Zhao,  “An  Adaptive  Control  for  Wind 
Turbine,”  presented,  at  the  International  Conference  on  Renewable  Energies  and 
Power  Quality  05,  Zaragoza,  Spain,  March  2005. 


77 


[12]  R.  Pena,  J.  C.  Clare,  and  G.  M.  Asher,  “Doubly  Fed  Induction  Generator  Using 
Back-to-Back  PWM  Converters  and  its  Application  to  Variable  Speed  Wind- 
energy  Generation,”  IEE  Proc.-Electr.  Power  Appl.,  vol.  143,  no.  3,  pp.  231-241, 
May  1996. 

[13]  S.  K  Salman  and  B.  Badrzadeh,  “New  Approach  for  Modeling  Doubly-Fed 
Induction  Generator  (DFIG)  for  Grid-Connection  Studies,”  presented  at  2004 
European  Wind  Energy  and  Exhibition,  London,  November  2004. 

[14]  F.  D.  Bianchi,  H.  De  Battista,  and  R.  J.  Mantz,  Wind  Turbine  Control  Systems,  pp. 
2-28,  Springer,  London,  2010. 

[15]  A.  I.  Munteanu,  A.  I.  Bratcu,  N-A  Cutululis,  and  E.  Ceang,  Optimal  Control  of 
Wind  Energy  Systems  Towards  a  Global  Approach,  Springer,  London,  2008. 

[16]  M.  O.  L.  Hansen,  Aerodynamics  of  Wind  Turbines,  Second  Edition,  Earthscan, 
London,  2008. 

[17]  A.  Betz,  Introduction  to  the  Theory  of  Flow  Machines,  Permagon  Press,  Oxford, 
1966. 

[18]  J.  F.  Manwell,  J.  G.  McGowan,  and  A.  L.  Rogers,  Wind  Energy  Explained: 
Theory,  Design,  and  Application,  John  Wiley  and  Sons  Ltd.,  West  Sussex, 
England,  2002. 

[19]  P.  C.  Krause,  O.  Wasynczuk,  and  S.  D.  Sudhoff,  Analysis  of  Electric  Machinery 
and  Drive  Systems,  Second  Edition,  John  Wiley  &  Sons,  2002. 

[20]  R.  De  Doncker,  D.  W.  J.  Pulle,  and  A.  Veltman,  Advanced  Electrical  Drives, 
Analysis,  Modelling,  Control,  Springer,  London,  2011. 

[21]  A.  L.  Julian,  Notes  for  EC4130  (Advanced  Electrical  Machinery  Systems),  NPS, 
2010  (unpublished). 

[22]  P.  Ioannou  and  J.  Sun,  Robust  Adaptive  Control,  Prentice  Hall,  1996. 

[23]  K.  J.  Astrom  and  B.  Wittenmark,  Adaptive  Control,  Second  Edition,  Dover 
Publications,  New  York,  2008. 

[24]  S.  Sastry  and  M.  Bodson,  Adaptive  Control,  Stability,  Convergence,  and 
Robustness,  Prentice  Hall,  New  Jersey,  1989. 

[25]  J.  J.  E.  Slotline  and  W.  Li,  Applied  Nonlinear  Control,  Prentice  Hall,  New  Jersey, 
1991. 


78 


[26]  K.  J.  Astrom  and  B.  Wittenmark,  Computer  Controlled  Systems,  Theory  and 
Design ,  Third  Edition,  Prentice  Hall,  New  Jersey,  1997. 

[27]  R.  G.  Brown  and  P.  Y.C.  Hwang,  Introduction  to  Random  Signals  and  Applied 
Kalman  Filtering,  Third  Edition,  Wiley,  New  York,  1997. 

[28]  J.  G.  Massey,  “Doubly  Fed  Induction  Machine  Control  for  Wind  Energy 
Conservation,”  M.S.  Thesis,  NPS,  Monterey,  CA,  2009. 

[29]  J.  H.  Yang,  “Robust  Adaptive  Control  Using  a  Filtering  Action,”  PhD 
dissertation,  NPS,  Monterey,  CA,  2009. 

[30]  C.  Rohrs,  L.  Valavani,  and  M.  Athans,  and  G.  Stein,  “Robustness  of  Continuous 
Time  Adaptive  Control  Algorithms  in  the  Presence  of  Unmodeled  Dynamics,” 
IEEE  Transactions  on  Automatic  Control,  vol.  30,  no.  9,  pp.  881-889,  September 
1985. 

[31]  N.  Hovakimyan  and  C.  Cao,  LI  Adaptive  Control  Theory,  Society  for  Industrial 
and  Applied  Mathematics,  Philadelphia,  2010. 

[32]  National  Oceanic  and  Atmospheric  Administration,  National  Weather  Service, 
Definition  of  Wind  Gust.  [Online].  Available: 
http://www.weather.eom/glossary/g.html  (accessed  August  9,  2011). 

[33]  R.  Cristi,  Notes  for  EC4300  (Advanced  Topics  in  Control  Systems:  Wind  Turbine 
Control),  NPS,  2010  (unpublished). 


79 


THIS  PAGE  INTENTIONALLY  LEFT  BLANK 


80 


INITIAL  DISTRIBUTION  LIST 


1 .  Defense  Technical  Infonnation  Center 
Ft.  Belvoir,  Virginia 

2.  Dudley  Knox  Library 
Naval  Postgraduate  School 
Monterey,  California 

3.  Dr.  Clark  Robertson 

Electrical  Engineering  and  Computer  Department 
Code  EC/R1 

Naval  Postgraduate  School 
Monterey,  California 

4.  Dr.  Roberto  Cristi 
Code  EC/Cl 

Naval  Postgraduate  School 
Monterey,  California 

5.  Dr.  Alexander  Julian 

Electrical  Engineering  and  Computer  Department 

Code  EC/J1 

Naval  Postgraduate  School 
Monterey,  California 

6.  Dr.  Xiaoping  Yun 

Electrical  Engineering  and  Computer  Department 

Code  EC/Yx 

Naval  Postgraduate  School 
Monterey,  California 

7.  Embassy  of  Greece 
Office  of  Naval  Attache 
Washington,  District  of  Columbia 

8.  LT  Stamatios  Orfanos  -  Pepainas 

Hellenic  Navy  General  Staff 
Athens,  Greece 


81 


